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A  BASIC  APPROACH  TO  THE  EVALUATION  OF 
RISK?  INTERRELATED  INVESTMENTS 
by 

Frederick  S.  Hillier 


1.  Introduction 

When  evaluating  proposals  for  capital  investment,  it  often  is 
necessary  to  consider  interrelationsh. ps  of  various  kinds  between  the 
projects.  For  example,  some  projects  may  be  complementary,  such  as 
those  which  would  benefit  in  common  from  the  facilities  that  would 
need  to  be  provided  for  any  one  of  them.  Some  proposals  may  even  be 
necessarily  contingent  on  undertaking  one  or  more  prerequisite  projects 
or  on  the  occurrence  of  particular  favorable  conditions.  Certain 
proposals  may  instead  be  mutually  exclusive  in  that  they  involve 
different  ways  of  performing  the  same  service.  Other  proposed  invest¬ 
ments  may  at  least  be  competitive  in  the  sense  that  undertaking  one 
reduces  the  benefits  that  would  result  from  also  undertaking  any  of 
the  others.  Furthermore,  some  or  all  of  the  proposals  may  need  to 
compete  for  certain  limited  resources  such  as  capital,  manpower,  and 
facilities. 

Fortunately,  considerable  progress  has  been  made  in  recent  years 
in  developing  ways  of  taking  these  interrelationships  into  account. 

Of  particular  note  is  Weingartner 1 s  comprehensive  work  [17]  on  applying 
mathematical  programming  to  the  analysis  of  such  capital  budgeting 
problems.-^  All  of  the  interrelationships  discussed  above  usually  can 

—  Also  see  the  survey  paper  by  Weingartner  [ 16 ] . 
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be  formulated  quite  conveniently  in  either  a  linear  or  nonlinear  pro¬ 
gramming  format,  so  this  provides  a  powerful  tool  for  the  analysis  of 
interrelated  projects. 

Another  crucial  factor  to  consider  when  evaluating  proposed  invest¬ 
ments  is  the  element  of  risk.  Unfortunately,  risk  is  a  relatively 
difficult  factor  to  deal  with  in  an  explicit,  quantitative  fashion. 
However,  Hillier  [8]  has  introduced  the  concept  of  analyzing  the  various 
individual  determinants  of  project  worth  and  estimating  the  most  likely 
outcome  and  degree  of  uncertainty  associated  with  each  in  order  to 
obtain  the  overall  probability  distribution  of  a  measure  of  merit  for 
the  investment  (e.g.,  distribution  of  internal  rate  of  return),  and 

then  using  this  distribution  as  a  basis  for  management's  evaluation 
2  / 

of  risk.—'  He  also  presented  an  analytical  approach  for  developing 
this  distribution.  Hertz  [k]  then  pointed  out  that  computer  simulation 
also  could  be  used  to  develop  this  same  information.  This  approach 
to  risk  analysis  has  been  adopted  quite  widely  during  the  last  few 
years. 

Thus,  by  using  the  approaches  mentioned  above,  it  is  now  possible 
to  deal  separately  with  (l)  the  effect  of  interrelationships  between 
proposed  projects,  and  (2)  tne  evaluation  of  the  risk  associated  with 
individual  investments.  However,  this  does  not  answer  the  question  of 
how  to  simultaneously  take  both  considerations  into  account  in  order  to 
choose  the  best  overall  combination  from  the  group  of  proposed  invest¬ 
ments.  Is  it  possible  to  combine  the  two  approaches  in  some  reasonable 
way  to  accomplish  this? 

-  Also  see  [7]  and  [10]. 
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Unfortunately,  this  kind  of  extension  usually  would  not  be  parti¬ 
cularly  feasible  with  the  simulation  approach  to  risk  analysis.  As 
the  author  points  out  elsewhere  [$,  pp.  84f. ],  this  approach  "has  a 
number  of  serious  disadvantages  which  make  it  poorly  suited  for  the 
analysis  of  risky  interrelated  investments.  One  of  the  lesser  of 
these  is  that  simulation  is  inherently  an  imprecise  technique,  even 
with  respect  to  the  model  used,  since  it  provides  only  statistical 
estimates  rather  than  exact  results.  In  other  applications,  these 
estimates  commonly  are  of  the  mean  of  a  distribution,  in  which  case 
the  ones  yielded  by  the  usual  computer  run  sizes  tend  to  have  a  low 
but  tolerable  precision.  However,  estimates  of  probability  distribu¬ 
tions  are  considerably  more  crude,  especially  wh^n  a  tail  of  the  dis¬ 
tribution  is  critical,  as  is  the  case  here.  Estimates  of  differences 
between  alternatives  also  are  at  least  as  imprecise.  Furthermore, 
simulation  is  a  cumbersome  way  to  study  a  problem,  since  it  requires 
developing  the  model  and  input  data,  and  then  doing  the  computer  pro¬ 
gramming  and.  executing  the  computer  runs.  However,  perhaps  the  most 
critical  disadvantage  is  that  simulation  yields  only  numerical  data 
about  the  predicted  performance  of  investments,  so  that  it  yields  no 
additional  insight  into  cause-and-effect  relationships.  Therefore, 
every  slightly  new  case  must  be  completely  rerun.  This,  in  addition 
to  the  imprecision,  tends  to  make  it  impossible  to  conducx  a  satis¬ 
factory  sensitivity  analysis.  An  even  more  serious  implication  is 
that  a  new  simulation  run  is  required  for  each  combination  of  inter¬ 
related  investments  under  consideration.  Since  even  one  run  of 
acceptable  length  is  expensive,  and  there  may  be  thousands  or  even 


3 


I 


millions  of  feasible  combinations,  this  too  would  tend  to  be  pro¬ 
hibitively  expensive.  Finally,  even  if  one  were  able  to  obtain  a 
crude  estimate  of  the  probability  distribution  of  present  value  for 
each  of  these  feasible  combinations,  the  simulation  approach  provides 
no  guidance  on  how  to  use  all  of  these  masses  of  data  in  order  to 
select  the  investments  to  be  approved." 

On  the  other  hand,  the  analytical  approach  to  risk  analysis  is 
relatively  well-suited  for  extension  to  the  case  of  interrelated 
investments.  In  fact,  it  can  be  incorporated  directly  into  a  mathe¬ 
matical  programming  format.  One  way  of  doing  this  is  by  means  of  a 
chance-constrained  programming  formulation.  This  method  has  been 
explored  in  some  depth  by  Byrnes,  Charnes,  Cooper,  and  Kortanek  [1,2], 
Wftslund  [13,l4],  Hillier  [5i  9t  Ch.  73 >  and  others.  An  alternative 
approach  which  may  provide  a  more  fundamental  and  precise  evaluation 
of  the  overall  merit  of  the  set  of  approved  investments  also  has  been 
developed  by  Hillier  [9]  in  a  companion  volume.  This  method  uses 
present  value  (treated  as  a  random  variable)  and  expected  utility  of 
present  value  as  criteria  in  order  to  choose  the  best  overall  combi¬ 
nation  from  the  group  of  proposed  investments. 

In  this  companion  volume  [ 9 ] j  the  investigation  begins  by  deriving 
properties  that  characterize  an  optimal  combination  of  investments.  A 
model  is  then  developed  for  describing  the  interrelated  cash  flows 
generated  by  a  given  set  of  investments.  This  model  would  be  used 
for  determining  the  mean,  variance,  and  possibly  the  functional  form 
of  the  probability  distribution  of  present  value.  The  question  of 
when  this  distribution  would  be  normal  or  approximately  normal  (even 


when  the  distributions  of  the  individual  cash  flows  are  not)  is  explored 
in  detail.  Next,  several  convenient  models  of  utility  functions  having 
desirable  properties  are  formulated.  For  each,  an  expression  is  derived 
for  calculating  the  expected  utility  of  approving  any  particular  combi¬ 
nation  of  investments.  Given  these  results,  an  approximate  linear 
programming  approach  and  an  exact  branch-and-bound  algorithm  are 
developed  for  determining  the  best  combination.  Numerous  suggestions 
regarding  the  practical  implementation  of  the  theory  and  procedures 
also  are  given. 

The  purpose  of  this  paper  is  to  focus  on  a  particular  basic  model 

ftrom  [9]-/  that  seems  to  be  particularly  well-suited  for  practical  use, 

and  to  elaborate  further  on  how  to  formulate  and  apply  it.  However, 

rather  than  repeating  most  of .the  technical  details  involved  in 

developing  this  model  and  its  supporting  theory,  the  emphasis  here  is 

on  motivating,  identifying,  and  interpreting  the  basic  structure  of 

the  model.  Thus,  this  is  to  be  a  relatively  nontechnical  expository 
4/ 

treatment,- 

In  pursuit  of  this  objective,  the  next  section  develops  the 
portion  of  the  r  odel  designed  for  considering  interrelationships 
between  investments.  Section  5  then  introduces  the  evaluation  of  risk 
into  this  framework.  Section  4  reviews  solution  procedure*  for  this 
model  that  were  developed  in  [9]  for  seeking  the  best  combination  of 
investments.  Computational  experience  with  these  procedures  is 
outlined  in  Section  5>  with  conclusions  given  in  Section  6. 

X] - - - 

~  See  primarily  Sections  1.6,  5-2,  5-3,  5-^,  and  6.1. 

4  / 

~  For  ease  of  exposition,  the  problem  is  discussed  in  terms  of  a 
business  firm  which  currently  has  a  number  of  opportunities  for  major 
capital  investments  or  projects  to  be  considered  by  top  management. 
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2.  A  Model  for  Considering  Interrelationships 

Suppose  that  a  number  of  proposals  for  capital  investment  have 
been  submitted  to  management  for  their  consideration.  As  discussed 
above,  it  is  likely  that  some  of  the  proposed  projects  are  interrelated 
in  one  or  more  significant  ways.  Therefore,  it  is  important  to  recog¬ 
nize  the  nature  of  these  interrelationships  and  to  take  their  effects 
into  account  in  the  analysis.  A  model  for  doing  this  is  developed 
below. 

The  first  assumption  is  that  all  decisions  to  be  made  are  of  the 
"yes-or-no"  rather  than  the  "how  much"  type.  Thus,  the  kind  of  invest¬ 
ments  being  considered  are  capital  projects  to  be  approved  or  rejected 

rather  than,  for  example,  common  stock  where  the  decisions  are  how 

5  / 

many  shares  of  each  type  to  purchase.-  However,  it  should  be  empha¬ 
sized  that  "investment  decision"  may  be  defined  in  a  broad  sense  here 
so  as  to  include  various  investment  strategy  possibilities.  Thus,  in 
addition  to  a  flat  approval  or  rejection  of  a  project,  other  alterna¬ 
tives  such  as  "postpone  the  investment  for  a  year"  or  "try  a  pilot 
run  first"  may  also  be  considered  by  treating  all  of  the  possibilities 
as  mutually  exclusive  investments  to  be  approved  or  rejected.  Similarly, 
several  alternative  levels  at  which  to  undertake  a  project  may  be 
treated  as  mutually  exclusive  investments.  Furthermore,  alternative 
strategies,  such  as  "try  project  A  for  a  year  and  then  switch  to 
project  B  if  the  losses  are  more  than  x  dollars"  for  several  different 
values  of  x,  also  are  mutually  exclusive  investments  which  may  be 

- — 

-  The  latter  problem  of  portfolio  selection  has  been  analyzed  effec¬ 
tively  by  Markowitz  [12]  and  others. 
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contingent  upon  certain  other  decisions  (such  as  not  undertaking  project 
B  immediately).  As  these  examples  illustrate,  confining  the  analysis 
to  yes-or-not  decisions  need  not  be  particularly  restrictive  in  a 
capital  budgeting  context. 

Now  consider  how  to  begin  formalizing  the  problem  mathematically. 

Suppose  that  all  of  the  interesting  possible  investment  decisions,  in 

the  broad  sense  described  above,  have  been  identified.  Let  the  number 

of  such  investment  decisions  be  denoted  by  m.  Then  introduce  a 

decision  variable  b,  for  each  of  the  decisions,  where  the  variable 

k 

is  assigned  a  value  of  one  or  zero  according  to  whether  a  decision  is 

"til 

if  the  k  investment  is  approved 
th 

if  the  k  investment  is  rejected  , 

for  k  •-  1,2,. .  ,,m.  Let  6  =  (b-^bg, . . .  ,6m) .  Hence,  a  "solution"  to 
this  capital  budgeting  problem  corresponds  to  a  particular  value  of  b 
where  each  of  the  components  is  either  zero  or  one.  However,  not  all 
combinations  of  zeroes  and  ones  need  be  a  "feasible"  solution.  Thus, 
if  a  particular  group  of  investments  are  mutually  exclusive,  then  the 
constraint , 


yes  or  no.  Thus, 


must  be  satisfied,  tf'nere  the  summation  is  over  the  indices  corresponding 
to  the  set  K  of  these  mutually  exclusive  alternatives.  One  such 
constraint  would  need  to  be  imposed  for  each  set  of  mutually  exclusive 
investments.  Any  other  restrictions  also  must  be  taken  into  account. 


7 


For  example,  if  investment  j  is  contingent  upon  investment  i  being 
approved,  this  can  be  expressed  by  the  constraint, 

6.  <  6.  , 

J  -  i 

since  then  8^  can  equal  one  only  if  6^  is  equal  to  one.  It  is  also 

common  to  have  budget  restrictions  whereby  the  total  outlay  for  all 

approved  projects  during  a  certai;,  period  or  periods  must  fall  within 

prescribed  limits.  Similarly,  limitations  on  required  resources  such 

as  labor  and  materials  may  rule  out  certain  combinations  of  investments. 

These  and  similar  restrictions  usually  can  be  expressed  conveniently 

in  a  mathematical  programming  format,  as  described  in  some  detail  by 

Weingartner  [17,  Sections  2.3,  3-6,  3*7>  7- 1 >  7-2]. 

In  order  to  systematically  determine  which  particular  feasible 

solution  8  is  optimal  (with  respect  to  the  model  being  developed), 

it  is  necessary  to  quantify  the  effect  of  a  choice  of  5.  In  very 

general  terms,  the  basic  effect  of  approving  a  particular  combination 

of  investments  is  to  generate  a  stream  of  cash  flows  or  imputed  cash 

flows.  Immediate  cash  outlays  (negative  cash  flows)  normally  are 

required  for  the  investments,  and  such  outlays  may  need  to  be  continued 

for  some  period  of  time..  The  pa,’ off  then  comes  in  the  form  of  income 

(positive  cash  flow)  over  some  interval  or  instant  of  time,  or  in  the 

form  of  other  benefits  of  real  value  to  the  investors.  In  order  to 

have  a  common  measure  for  these  consequences,  the  other  benefits  will 

be  measured  in  terms  of  their  equivalent  cash  value  (positive  imputed 

,  6/ 

cash  flow).- 

— ■  Hereafter,  the  term  "cash  flow"  will  be  taken  to  encompass  both  real 
cash  flow  and  imputed  cash  flow.; 


For  purposes  of  analysis,  the  future  is  divided  into  time  periods 
of  convenient  length  (e.g.,  a  year).  Let  n  be  the  maximum  number  of 


periods  into  the  future  in  which  significant  cash  flows  may  occur 
because  of  the  investments.  Immediate  cash  flows  are  considered  to 
take  place  in  period  0.  Thus,  it  is  necessary  to  analyze  what  the 
cash  flows  will  be  in  periods  0,l,2,...,n.  Let  X.(8)  be  the  total 
net  cash  flow  (total  positive  cash  flow  minus  total  negative  cash  flow) 
during  time  period  j  that  would  result  from  making  the  particular 
decision  6  (j  =  0,1, ...,n). 

It  is  important  to  recognize  two  facts  about  the  nature  of  the 

X.(8).  The  first  is  that  the  cash  flow  stream  resulting  from  a  decision 
J 

8  is  actually  an  aggregation  of  many  distinct  cash  flow  streams,  some 
of  which  are  interrelated  in  various  ways.  Each  individual  investment 
approved  generates  a  cash  flow  stream,  which  may  itself  be  an  aggre¬ 
gation  of  cash  flow  streams  from  a  number  of  sources  representing  the 
different  kinds  of  outlays  or  incomes  for  the  investment.  However,  the 
numerical  values  associated  with  the  individual  cash  flow  streams  may 
be  affected  significantly  by  the  others.  It  is  the  analysis  of  these 
interrelationships  between  the  components  of  the  aggregated  cash  flow 
stream  that  forms  the  basis  for  the  model  developed  in  this  section. 

The  second  crucial  property  of  the  X.(8)  is  that,  when  dealing 

J 

with  risky  investments,  they  normally  are  random  variables  rather  than 
constants  (with  the  probable  exception  of  Xq(8)).  It  is  inherent  in 
the  nature  of  risky  investments  that  there  is  considerable  uncertainty 
about  their  outcomes.  Thus,  a  decision  8  can  result  in  the  total 
net  cash  flow  during  period  j  being  any  one  of  many  different  possible 


values,  depending  on  intervening  circumstances.  Therefore,  the  value 
that  X.(6)  will  actually  take  on  can  only  be  described  in  terms  of 
a  probability  distribution. 

Thus,  the  direct  impact  of  a  decision  5  is  that  an  aggregate 

cash  flow  stream  is  generated  from  the  joint  probability  distribution 

of  Xq(6),X1(5) , . . . ,Xn(s) .  In  order  to  evaluate  a  particular  5, 

and  to  choose  between  alternative  feasible  values  of  6,  it  is 

therefore  necessary  to  evaluate  the  desirability  of  different  cash 

flow  streams.  Although  various  measures  of  desirability  have  been 

used  in  practice,  there  is  considerable  theoretical  support  for  using 

discounted  cash  flow  methods  and  particularly  the  present  value  method. 

The  latter  criterion  is  the  one  adopted  here.  Thus,  let  i  be  the 

rate  of  interest,  normally  the  market  cost  of  capital,  which  properly 

7/ 

reflects  the  invesuor's  time  value  of  money.-  Then  let  P(5)  be 
the  total  net  present  value  of  the  investments  approved  by  5,  so 
that 


P(5) 


n  X.(6) 

2?  7^- 

j«o  (1+1 r 


Since  the  X.(5)  are  net  cash  flows  (income  minus  outlays)  for  the 
0 

respective  periods,  P(8)  thereby  represents  the  profit  (if  positive) 
or  loss  (if  negative)  from  the  approved  investments  after  discounting 
for  the  time  value  of  the  money  invested. 

I-,  should  be  emphasized  that,  since  the  X.(5)  are  random 
variables,  P(5)  also  is  a  random  variable.  When  a  set  of  risky 

—  This  interest  rate  may  be  set  at  different  values  in  different 
periods  if  desired. 
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investments  is  approved,  it  is  not  definitely  known  in  advance  whether 
they  will  result  in  a  large  loss  or  a  large  profit  or  something  in 
between.  Therefore,  "the"  present  value  P(6)  that  will  result  from 
a  decision  6  can  only  be  represented  in  terms  of  a  probability  dis¬ 
tribution  over  its  range  of  possible  values,  as  shown  in  Figure  1. 

When  only  a  single  investment  is  being  evaluated,  the  essence  of 
the  risk  analysis  approach  proposed  by  Hillier  [8]  and  Hertz  [4]  is 
to  develop  such  a  probability  distribution  (for  either  present  value  or 
another  measure  of  merit)  and  then  to  analyze  this  "risk  profile"  in 
order  to  conclude  whether  the  amount  of  risk  is  acceptable  or  not. 

If  the  degree  of  risk  represented  by  the  probabilities  and  magnitudes 
of  possible  losses  is  sufficiently  small  relative  to  the  values  of 
these  quantities  for  the  possible  gains,  then  the  investment  would  be 
approved.  For  the  problem  now  under  consideration,  where  many  inter¬ 
related  investments  need  to  be  evaluated  simultaneously,  the  purpose 
of  the  analysis  changes  somewhat.  Rather  than  assessing  whether  a 
particular  risk  profile  is  "acceptable,"  the  main  objective  becomes 
to  identify  the  "best"  obtainable  risk  profile.  In  effect,  the  goal 
is  to  determine  which  combination  of  investments  provides  the  "best" 
probability  distribution  of  P(6),  considering  both  the  possibilities 
of  losses  and  of  gains.  An  approach  to  making  this  kind  of  choice 
between  alternative  probability  distributions  is  developed  in  the 
next  section.  Meanwhile,  the  question  remains  of  how  to  analytically 
determine  the  probability  distribution  of  P(5)  for  any  given  value 
of  5,  taking  into  account  the  role  played  by  the  various  kinds  of 
interrelationships  between  the  approved  investments.  This  is  the 
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question  that  will  be  explored  throughout  the  remainder  of  this 
section. 

Before  considering  interrelationships  between  investments,  it  is 
necessary  to  investigate  the  individual  investment  proposals  by  them¬ 
selves.  Thus,  the  first  step  in  determining  the  distribution  of  P(8) 
is  to  consider  each  investment  in  isolation,  analyzing  what  its 
performance  would  be  if  it  were  the  only  one  approved  (except  for 
any  necessary  predecessors).  To  illustrate,  consider  a  typical  invest¬ 
ment  k  (k  r-  l,2,...,m)  in  isolation.  This  investment  generates  one 
or  more  cash  flow  streams  from  its  sources  of  outlays  or  incomes. 

Denote  the  present  value  of  the  aggregation  of  these  cash  flow  streams 
(when  no  other  investments  are  made)  by  P^.  If  this  investment  has 
any  risk  associated  with  it,  then  P^  actually  is  a  random  variable 
having  some  underlying  probability  distribution.  Therefore,  the  cash 

flow  streams  should  be  analyzed  in  order  to  estimate  at  least  the  mean 

2  8/ 

and  variance  a of  this  distribution.-  The  expected  present 

value  nk  is  merely  the  sum  of  all  of  the  expected  discounted  cash 

flows  over  the  respective  time  periods  (0,1,..., n).  Determining  the 
2 

variance  cr^  is  somewhat  more  difficult,  particularly  because  some 
the  cash  flows  may  be  correlated.  One  fairly  likely  kind  of  corre¬ 
lation  is  between  cash  flows  from  the  same  source  in  different  periods. 
For  example,  it  is  quite  common  that  if  an  investment  performs  con¬ 
siderably  better  (or  worse)  than  expected  in  the  early  periods,  then 
it  will  tend  to  also  perform  better  (or  worse)  in  subsequent  periods 

- - 

-  The  estimation  of  these  and  other  parameters  of  the  model  is  dis¬ 
cussed  at  the  end  of  the  section,  and  references  giving  practical 
estimating  procedures  are  cited. 
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than  had.  been  expected  initially.  The  degree  of  this  tendency  can  be 
measured  by  the  correlation  coefficient  between  the  cash  flows  in  at 

least  consecutive  time  periods  (as  described  further  in  [9,  Sec.  A.2]). 

2 

Some  models  for  calculating  in  such  cases  are  described  in  [8] 
and  illustrated  in  [lOj.  The  second  (and  perhaps  less  likely)  kind  of 
correlation  is  between  different  cash  flow  streams.  The  approach  in 
this  case  is  completely  analogous  to  that  described  subsequently  for 
considering  the  correlation  between  the  aggregate  cash  flow  streams  of 
different  investments  in  order  to  calculate  the  standard  deviation  of 
P(5).-// 

2 

Suppose  that  one  has  estimated  the  mean  ^  and  variance  a  of 

K  k 

the  present  value  of  investment  k,  assuming  it  is  the  only 

investment  approved,  for  each  k  =  l,2,...,m.  If  there  were  no  inter¬ 
relationships  of  any  kind  between  the  investments,  it  would  then  be 
necessary  only  to  add  these  values  for  the  approved  investments  in 
order  to  find  the  mean  and  variance  of  P(6)„  However,  the  presence 
of  certain  kinds  of  interrelationships  can  invalidate  this  simple 
additive  relationship.  Consider  first  how  this  can  happen  with  the 
mean  because  of  complementary  or  competitive  effects.  Two  investments 
may  be  complementary  because  they  would  share  common  costs  or  be 
mutually  reinforcing  in  generating  income.  Therefore,,-  even  if  either 
investment  by  itself  would  be.  unattractive,  the  combination  of  the  two 
investments  together  might  be  very  worthwhile.  Conversely,  two  com¬ 
petitive  investments  might  be  very  attractive  on  an  individual  basis,  but 
still  bf  very  undesirable  m  combination  In  both  cases,  what  is  happening, 
in  effect,,  is'  that  t h<  pr..,mt  valu-  of  both  investments  together  is 


something  different  than  the  sum  of  the  present  values  that  each  would 
attain  if  the  other  were  not  approved.  Thus, 


m 

P(8)  =  £  P.  8  +  h(5) 
k=l  K  K 

where  the  function  h(5)  is  the  net  amount  by  which  P(&)  needs  to 
be  adjusted  due  to  complementary  or  competitive  interactions  between 
the  approved  investments.  In  general,  it  is  possible  to  have  the 
amount  of  a  complementary  or  competitive  effect  be  a  random  variable 
rather  than  a  constant,  and  to  have  some  of  the  interactions  involve 
more  than  two  investments  simultaneously  in  a  complicated  way.  However, 
for  simplicity,  the  model  here  assumes  that  all  these  effects  are  both 
deterministic  and  "pairwise  additive,"  so  that  any  joint  effect 
involving  more  than  two  investments  is  merely  the  cumulation  of  the 
pairwise  effects.  As  a  result. 


h(6)  = 


m 


m 


2  2 
j=l  k=l 
k^j 


VA 


where  (u,,  +  u..)  is  the  net  addition  (positive  or  negative)  to  total 
jk  kj 

present  value  due  to  complementary  or  competitive  interactions  (if  any) 

between  investments  j  and  k  if  both  are  made.  Although  only  the 

sum  (u,,  +  u,  ,)  is  relevant,  the  convention  is  adopted  here  that 
jk  kj 

u  =  u  , ,  so  Lij,  represents  investment  j's  "equal  share"  of  this 
jk  kj'  jk 

effect.  Thus,  whereas  |ij  or  would  be  the  individual  contribu¬ 

tion  to  expected  total  present  value  of  investment  j  or  k  made  by 
itself,  their  joint  contributions  become  (mj  +  U^)  and  (nk  + 
if  both  are  approved.  Each  investment  j  may  have  a  nonzero 
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with  more  than  one 


complementary  or  competitive  interaction 
investment  k.  Therefore,  letting  p(6) 


denote  the  expected  value 


of  P(8),  this  quantity  becomes  simply 


Now  consider  how  to  find  the  varianc  of  P(5),  which  will  be 
2 

denoted  by  a  (6).  Since  h(6)  is  assumed  to  be  deterministic  in  the 
expression  for  P(5)  given  above,  complementary  and  competitive,  inter¬ 
actions  have  no  effect  on  this  variance.-  However,  there  may  well  be 
other  kinds  of  interrelationships  between  the  investments  that  affect 
the  amount  by  which  the  realized  value  of  P(5)  deviates  from  its 
expectation.  In  particular,  any  common  factors  that  help  determine 
the  actual  performance  of  the  investments  relative  to  expectations 
are  of  this  kind.  These  factors  may  be  internal  to  the  firm,  such  as 
the  outcome  of  labor  negotiations  or  the  time  at  which  resources  from 
an  ongoing  project  will  become  available  for  reallocation.  However, 
the  most  important  factors  probably  are  exogenous,  such  as  the  general 
state  of  the  economy  or  advances  made  by  competitors.-  If  there  are 
any  such  factors  that  may  exert  a  general  influ-..ce  on  the  performance 

of  some  of  the  investments,  the  consequence  is  chat  the  P.  (j  = 

J 

l,f,...,m)  would  be  correlated  rather  than  statistically  independent. 


Let  p.,  denote  the  correlation  between  P.  and  P,  ,  and  let 
Jk  ,1  k' 

a  =  p ,,oo,  be  the  corresponding  covariance.  Therefore,  assuming 
jk  jk  j  x 

that  estimates  of  these  parameters  can  be  obtained  (as  will  be  dis¬ 


cussed  shortly),  a  (6)  can  be  calculated  simply  as 


1L 


o2(  6) 


+ 


m 


S 


o..  6. 

jk  k 


5.  . 

J 


2 

Given  the  mean  p(5)  and  variance  o  (8)  of  total  present  value 
P(&),  the  only  remaining  question  regarding  the  probability  distri¬ 
bution  of  this  random  variable  is  its  shape  (functional  form).  One 

simple  case  is  where  the  P.  are  normally  distributed  (i.e., 

0 

^l,^2,'"‘,^m  ^as  a  mu  bivariate  normal  distribution),  so  that  the 

sum  yielding  P(5)  also  has  a  normal  distribution.  Unfortunately, 

the  distribution  of  the  P.  often  would  differ  significantly  from 

J 

normal,  and  most  other  distributions  are  not  preserved  under  addition. 
The  nature  of  risky  investments  is  such  that  their  distributions  fre¬ 
quently  are  skewed  to  the  left,  as  depicted  in  Figure  1,  or  even 
bimodal.  In  such  cases,  it  is  not  possible  to  draw  any  conclusion 
as  to  the  exact  functional  form  of  the  distribution  of  P(&).  However, 
there  is  a  good  chance  that  an  approximate  answer  can  be  given.  To 
motivate  this,  recall  that  P(5)  is  the  discounted  aggregation  of 
many  different  cash  flow  streams,  with  perhaps  several  such  streams 
being  generated  by  each  of  the  approved  investments.  Then  recall  that 
the  Central  Limit  Tneorem  indicates  that,  when  summing  random  variables 
having  distributions  other  than  normal,  the  distribution  f  this  sum 
still  approaches  a  normal  distribution  under  certain  conditions.  The 
conditions  required  actually  are  not  very  stringent.  Although  the 
most  commonly  known  version  of  the  Central  Limit  Theorem  requires  that 
the  random  variables  be  independent  and  identically  distributed, 
various  other  versions  also  exist  where  one  or  both  of  these  assumptions 
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This  has  two  important 


r */ 

can  be  replaced  by  much  weaker  conditions 
implications  for  the  current  model.  First,  since  the  present  value  of 
each  individual  cash  flow  stream  is  the  (discounted)  sum  of  the  cash 
flow  random  variables  for  the  respective  periods,  its  distribution 
actually  may  be  much  closer  to  normal  than  would  have  been  expected 
from  the  perhaps  highly  skewed  or  bimodal  distributions  for  these 
random  variables.  Second,  except  for  the  constant  h(&)  term,  p(8) 
is  the  sum  of  the  present  value  random  variables  for  the  individual 
cash  flow  streams  generated  by  the  approved  investments.  Therefore, 
the  distribution  of  P(8)  may  itself  be  much  closer  to  normal  than 
would  have  been  expected  from  the  distributions  of  the  individual 
present  values.  For  these  two  reasons,  it  often  would  be  a  reasonable 
approximation  to  take  the  distribution  of  P(8)  to  be  normal,  even 
when  the  distributions  for  individual  cash  flows  are  far  from  normal. 

Therefore,  if  it  is  feasible  to  estimate  the  parameters  of  the 
2 

model  (the  p.,  p.,  ,  cj.,  and  o  ),  then  the  problem  of  how  to 
J  0K  0 

identify  the  probability  distribution  of  present  value  P(&)  for  any 
given  feasible  solution  8  is  at  least  partially  solved.  The  mean 
p(8)  and  variance  o  (8)  can  be  calculated  by  the  expressions  given 
above,  and  the  form  of  the  distribution  can  perhaps  be  concluded  to  be 
approximately  normal.  As  will  be  seen  later,  this  much  information 
about  the  distribution  is  adequate  for  the  required  analysis. 

Hence,  it  is  crucial  to  be  able  to  estimate  the  parameters  of 
the  model.  Actually,  this  is  not  as  imposing  a  task  from  a  technical 


w; 


These  versions  are  presented  and  discussed  in  detail  in  the 


companion  volume  [9,  Sec.  4.2].- 
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standpoint  as  it  might  appear.  The  companion  volume  [9*  App. ]  describes 
simple  estimating  procedures  whereby  personnel  without  technical  back¬ 
grounds  can  provide  the  needed  information,  which  can  then  be  converted 
into  the  desired  estimates.  The  basic  approach  is  patterned  after  PERT. 
Thus,  Section  A.l  of  [9]  proposes  obtaining  not  one  but  three  estimates 
for  each  cash  flow,  namely,  a  "most  likely"  estimate,  an  "optimistic" 
estimate,  and  a  "pessimistic"  estimate.  These  three  estimates  can  then 
be  converted  into  estimates  of  the  mean  and  variance  of  the  cash  flow. 
Section  A. 2  presents  a  useful  model  for  correlation  patterns  between 
cash  flows  in  different  time  periods  and  between  different  cash  flow 
streams.  Section  A.b  indicates  how  the  three-estimate  approach  can  be 

extended  in  order  to  construct  estimates  of  correlation  coefficients. 

2 

Combining  these  tools  leads  to  the  desired  estimates  of  the 
and  o  by  only  obtaining  other  readily  comprehendable  types  of 
estimates.  Section  A. 6  discusses  the  estimation  of  complementary  or 
competitive  effects  (for  a  more  complicated  model).  For  most  purposes, 
it  should  be  adequate  to  use  only  "most  likely"  estimates  in  order  to 
construct  the  desired  estimates  of  the 

Thus,  the  model  for  considering  interrelationships  presented  in 
this  section  allows  one  to  systematically  identify  (at  least  in  part) 
and  compare  the  probability  distribution  of  present  value  P(5)  for 
the  various  feasible  solutions  5.  The  next  question  that  needs  to 
be  addressed  is  how  to  assess  the  element  of  risk  in  making  these 
comparisons.  A  model  for  this  purpose  is  developed  in  the  next 


section. 


3-  A  Model  for  Considering  Risk 

Suppose  that  the  probability  distribution  of  present  value  P(6) 
has  been  identified  for  a  number  of  different  feasible  solutions  8. 

For  example,  two  such  distributions  are  shown  in  Figure  2,  one  corres¬ 
ponding  to  a  fairly  conservative  set  of  investments  and  the  other  to  a 
relatively  risky  set.  It  is  somevrhat  more  likely  that  the  latter  set 
will  yield  the  larger  gain,  but  there  is  also  a  significant  probability 
that  it  will  yield  a  large -loss.  Therefore,  it  is  necessary  for 
management  to  assess  these  risk  profiles  and  evaluate  the  consequences 
of  the  different  possible  outcomes  in  order  to  choose  between  two  such 
sets  of  investments. 

In  general-,  there  would  be  many  alternative  sets  of  investments 
to  be  considered.  In  fact,  with  m  individual  investment  proposals, 
there  can  be  as  many  as  2m  feasible  solutions  5,  each  with  its 
own ^distribution  of  present  value.  It  is  necessary  to  choose  between 
all  of  these  distributions  in  order  to  determine  the  best  combination 
of  investments.  There  usually  would  be  far  too  many  interesting  com¬ 
binations  for  it  to  be  feasible  for  management  to  evaluate  and  compare 
all  of  these  alternatives  explicitly.  Therefore,  what  needs  to  be 
done  is  to  formalize  management ' s  judgment  on  risk  tradeoffs,  expressing 
it  quantitatively  so  that  the  selection  procedure  may  be  carried  out 
systematically  on  an  electronic  computer.— ^  Fortunately,  utility 
theory  provides  a  relatively  satisfactory  way  of  doing  just  this.  The 

I 

approach  involves  developing  management's  utility  function  U(p),  as 


In  actuality,  it  will  be  management,  not  a  computer,  that  will 


make  the  final  decision.  However,  the  computer  is  needed  to  identify 
the  few  best  alternatives  that  should  be  explicitly  evaluated  and 


compared  by  management. 


20 


illustrated  in  Figure  3,  where  U(p)  is  the  utility  if  p  is  the 

12/ 

realized  present  value  of  the  approved  set  of  investments. — ■  A  basic 
interpretation  of  utility  is  that  it  measures  the  relative  desirability 
(if  positive)  or  undesirability  (if  negative)  of  the  outcome.  Thus, 
suppose  that  p  =  1  corresponds  to  a  present  value  of  $1  (and  that 
the  slope  of  U(p)  remains  equal  to  one  for  0  <  p  <  l).— ■  For  any 

other  value  of  p,  U(p)  then  converts  a  present  value  of  that  amount 
to  its  equivalent  worth  relative  to  the  first  dollar  of  present  value 
in  terms  of  one's  willingness  to  expend  effort  and  take  risks  to 
achieve  it  (or  to  avoid  it  if  p  is  negative) . 

This  concept  of  U(p)  can  be  made  more  concrete  by  comparing 
risk-taking  situations.  For  example,  suppose  that  the  choice  is 
between  a  very  safe  and  a  relatively  risky  set  of  investments.  The 
safe  investments  will  yield  a  positive  present  value  of  exactly  p^. 

The  second  set  of  investments  would  have  a  50-50  chance  (probability 
of  l/2)  of  only  breaking  even  (p  =  0)  or  of  yielding  a  present 
value  of  exactly  p2-  Then  a  preference  for  the  safe  set  of  invest¬ 
ments  implies  that  ^  U(p2)  <  U(p1),  whereas  a  preference  for  the 
second  set  implies  that  ^  U(pg)  >  U(p1).  Thus,  U(p2)  =  2U ( p1 )  is 
the  break-even  value  of  p2  (as  shown  in  Figure  3),  since  then  the 
expected  contribution  to  the  firm's  welfare  (i.e.,  expected  utility) 
would  be  exactly  the  same  for  the  two  alternatives. 

— 7  To  fix  the  scale  of  u(p),  the  convention  is  adopted  here  that 
this  function  passes  through  the  origin  with  slope  one. 

13/ 

—  For  ease  of  exposition,  money  often  will  be  assumed  here  to  be 
measured  in  units  of  dollars,  although  another  denomination  of  meaning¬ 
ful  size  or  a  denomination  from  another  monetary  system  could  of  course 
be  used  instead. 
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As  a  second  example,  suppose  that  the  choice  is  between  doing 
nothing  (so  p  =  0)  or  undertaking  a  risky  set  of  investments  that 
would  yield  either  a  negative  present  value  of  exactly  p  or  a 
positive  present  va^  .e  of  exactly  p+,  each  with  probability  -■ 

To  make  this  decision,  management  would  need  to  analyze  the  risk 
profile  of  the  set  of  investments,  and  decide  how  the  undesirability 
of  the  loss  p  compares  with  the  desirability  of  the  gain  p 

-  T- 

Approval  of  the  investments  implies  that  U(p+)  >  -U(p  ),  whereas 
rejection  implies  that  U(p+)  <  -U(p  ).  Thus,  the  value  of  p+  such 
that.  U(p+)  =  -U(p  )  is  the  break-even  value  (as  shown  in  Figure  3) 
above  which  the  risk  of  incurring  the  loss  p  is  justified  by  the 
largeness  of  this  prospective  gain. 

Thus,  the  utility  function  U(p)  merely  quantifies  management's 
judgment  on  risk  tradeoffs.  Given  this  function,  any  decision  6 
can  then  be  evaluated  by  calculating  its  expected  utility,  E{u(P(8))}, 
where  the  expectation  is  taken  with  respect  to  the  probability  dis¬ 
tribution  of  present  value  P(8).  This  expected  value  may  be  thought 
of  as  the  difference  of  two  quantities,  namely,  a  positive  term 
(expected  positive  utility)  built  up  from  the  gains  that  may  be 
obtained  minus  an  appropriate  penalty  for  the  risk  of  possible  losses. 
Hence,  E{u(P(8))]  is  a  valid  single-valued  measure  of  the  merit  of 
8  that  appropriately  combines  information  about  the  probable  gains 
and  the  risk  associated  with  the  approved  set  of  investments.  This 
reduces  the  problem  to  determining  which  feasible  combination  of 
investments  maximizes  expected  utility,  i.e., 

max  E{u(P(8) ) ) 


2k 


over  all  feasible  solutions  5.  Methods  for  doing  this  will  be 
discussed  in  the  next  section.  However,  before  turning  to  this  topic, 
the  crucial  question  of  how  to  develop  the  appropriate  utility  function 
U(p)  will  now  be  considered. 

Although  the  interpretation  of  u(p)  was  suggested  above,  it  is 

still  necessary  to  describe  just  what  information  is  needed,  how  to 

obtain  it  from  management  in  a  realistic  way,  and  how  to  use  it  to 

generate  U(p).  To  do  this,  two  alternative  models  for  the  shape  of 

the  utility  function  are  developed  below  which  require  obtaining  only 

three  relatively  nontechnical  decisions  from  management. 

How  should  utility  functions  U(p)  ordinarily  be  expected  to 

behave?  Several  characteristics  suggest  themselves  as  a  starting 

point  for  developing  a  model.  First,  U(p)  should  be  monotone 

increasing,  i.e.,  it  is  always  better  to  increase  present  value. 

Second,  since  tiny  changes  in  present  value  ordinarily  should  not 

matter  much,  it  is  likely  that  U(p)  would  be  continuous  and  possess 

continuous  derivatives  of  all  orders.  (The  n  derivative  will  be 

denoted  here  by  l/n^(p).)  Third,  it  is  plausible  that  U(p)  would 

be  a  concave  function,  since  the  incremental  worth  of  the  next  dollar 

of  present  value  should  not  increase  as  present  value  increases  (i.e., 

decreasing  marginal  utility).  Thus,  recalling  that  (0)  =  1,  the 

marginal  utility  U^(p)  would  always  be  decreasing  from  one  toward 

zero  as  p  increases  over  positive  values,  whereas  it  would  always 

14/ 

be  increasing  above  one  as  p  decreases  over  negative  values. — 

157 - 

— '  "decreasing"  and  "increasing"  are  used  here  in  the  weak  sense 
which  also  allows  "remaining  the  same." 
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Ihe  next  quest  lor  is  what,  happens  wnen  p  becomes  very  large  posi¬ 
tively  or  negatively.  For  the  positive  direction,  does  U^(p) 
continue  decreasing  all  the  way  to  zero  in  the  limit  or  does  it  instead 
converge  to  some  strictly  positive  lower  bound  as  p  increases? 
Whichever  the  case,  U(p)  thereby  converges  to  a  linear  asymptote 
ai  +  ^i?  (assuming  a^  is  finite),  where  b^  >  0  is  this  lower 
bound  on  l/^(p).  For  p  growing  very  large  in  the  negative  direc¬ 
tion  instead,  U^\p)  may  either  be  converging  to  some  upper  bound 
or  increasing  indefinitely  as  p  decreases.  In  the  former  case, 

U(p)  would  thereby  converge  to  a  linear  asymptote  a..,  -t-  b^P  (assuming 
is  finite)  as  p  -►  -  <*,  where  the  upper  bound  on  U^(p)  is 
b,,.  For  the  latter  possibility,  U(p)  would  decrease  with  p  at 
an  exponential  rate. 

Both  models  for  u(p)  have  all  of  the  properties  described 
above.  The  distinction  between  them  lies  in  which  possibility  occurs 
as  p  grows  very  large  in  the  negative  direction.  Thus,  the  first 
model  (which  will  be  called  the  "basic  model"  because  of  its  flexi¬ 
bility)  assumes  that  u(pj  converges  to  linear  asymptotes  in  both 
the  negative  and  positive  directions,  shown  in  Figure  4.  The 
well-known  class  of  functions  having  such  asymptotes  are  th<=  hyperbolas, 
and  it  is  further  assumed  that  u(p|  has  precisely  this  algebraic 
form.  The  parameters  of  the  asymptotes  arc-  assumed  to  be  finite,  and 

to  satisfy  the  conditions,  a  >  0,  0  <  b^  <  1,  a^  >  0,  and 

'l) 

b^  >  1.  Because  of  the  convention  that  u(0)  ~  0  and  Uv  (0)  -  1, 
these  parameters  are  further  required  i.o  satisfy  the  relation, 


u(p) 


Figure  4.  The  basic  model  for  U(p). 
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a  b  t  9  c  =  a  r  a  ,  ir.  order  for  the  desired  hyperbola  to  exist.— 7 
1  d.  c  1  1  C 

To  interpret  Tr.it  relit  ion,  let  d  aenote  the  number  such  that 
+  b^P  intersects  the  4-5°  line  through  the  origin  at  (-d,-d), 
as  shown  in  Figure  4,  Thus,  a.^  =  -d  ■+•  b^d  =  d(bp-l).  Since  the 

(l-b1)a2 

relation  can  also  be  written  as  a_  =  it  is  seen  that 

1  bg  -  ± 

a  =  d  ~  b,d,  so  that  a1  +  b  p  must  intersect  the  45°  lin>-  at 
L  i  i  1 

(dtU'  Therefore-  the  model  is,  in  effect,  defined  by  three  parameters- 

ct,  b- >  and  b^-vvie-r^  b^(0  <  b^  <  l)  is  the  slope  of  an  asymptote 

passing  through  (cl,d),  ^2^2  >  ^  bbe  sbope  °f  an  asymptote 

passing  through  C-d.-a),  and  U(p)  is  the  resulting  hyperbola. 

passing  through  the  origin  with  slope  one. 

The  expression  for  u(p)  defined  by  this  hyperbola  is  derived 

in  the  companion  volume  [9,  PP-  and  is  found  to  be 

(a  +  b  p]  +  (a  +  b  p)  -  Q 
U(p)  -  ^ ^ - 

where 

Q,  -Vi  (’a ^  bpP!  +  (a2  f  b2p)l?  "  ^ai  f  bib2P  +  a2^  * 

a1  =  d(l  -  b.)  , 
a 2  -  dfb^  -  1)  . 

In  the  next  section,  it  al_o  will  be  useful  to  know  the  first  thr.- 
derivatives  of  U(p). ,  which  are  easily  found  to  be 

— '  ibis  is  shown  in  the  companion  volume  [9>  P-  4l].; 
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where 


T  =  (b2  -  b1)2p  +  (a1  +  a2)(b1  +  bg  -  2)  . 

Thus,  in  order  to  construct  U(p)  for  the  basic  model,  it  is 
only  necessary  to  obtain  judgment  decisions  from  management  that  lead 
to  estimates  of  d,  b^,  and  bg.  The  question  to  be  considered  now 
is  how  to  interpret  each  of  these  three  parameters  of  the  model  in  a 
relatively  nontechnical  way,  so  that  the  needed  information  can  be 
obtained  from  management  in  terms  meaningful  to  them.  To  this  end, 
note  that  (-d,-d)  is  the  point  in  Figure  4  at  which  the  "piece-wise 
linear  approximation"  of  U(P),  (the  dashed-line  function  in  the 
figure),  changes  from  a  slope  of  b2  to  a  slope  of  one.  Hence, 
roughly  speaking,  the  additional  detrimental  effect  of  each  incremental 
dollar  (or  larger  denomination)  of  loss  is  comparable  to  that  for  the 
first  dollar  if  the  total  loss  is  somewhat  less  than  d,  whereas  it 
is  considerably  more  serious  if  the  total  loss  is  somewhat  larger  than 
d.  Thus,  when  examining  the  effect  of  increasingly  large  losses,  d 
is  the  magnitude  of  loss  above  which  the  situation  would  deteriorate 
at  an  accelerated  rate.  This  interpretation  of  d  can  be  explained 
to  management  in  terms  of  the  ability  to  absorb  additional  losses. 

Thus,  when  the  total  loss  (in  present  value)  is  relatively  small,  an 


29 


additional  dollar  (or  another  raeaningrully  large  denomination  of  money) 
of  loss  can  be  absorbed  almost  as  easily  as  the  first  one.  However, 
when  the  total  loss  is  already  relatively  large,  an  additional  loss 
would  be  much  more  serious  and  difficult  to  absorb.  What  is  the  point 
of  demarcation  between  these  two  situations?  How  large  can  the  loss 
be  before  it  has  become  much  more  difficult  to  absorb  any  additional 
loss?  After  a  suitable  explanation,  management's  answer  to  this  kind 
of  question  provides  the  desired  estimate  of  d. 

To  interpret  b^  and  b^,  recall  that  the  marginal  utility  (the 
slope  of  U(p))  is  essentially  one  when  p  is  close  to  zero,  essen¬ 
tially  b^  when  p  is  very  large,  and  essentially  b^  when  p  is 
a  very  large  negative  number.  Therefore,  b^  and  b^  essentially 
measure  the  relative  desirability  of  an  incremental  dollar  of  present 
value  when  a  very  large  gain  or  loss  had  already  been  incurred, 
measured  in  "equivalent  dollars"  of  worth  corresponding  to  the  desira¬ 
bility  of  an  incremental  dollar  when  the  realized  present  value  had 
been  close  to  zero.  To  make  this  more  concrete,  it  is  suggested  that 

X 

the  following  hypothetical  situation  be  posed  to  management.  Suppose 
that  a  project  has  already  been  approved  (to  the  exclusion  of  all  other 
proposed  investments)  that  will  lead  to  one  of  two  possible  outcomes. 
The  first  possibility  is  that  it  will  "break  even,"  yielding  a  present 
value  of  zero.  The  second  possible  outcome  is  that  the  project  will 
"make  it  big,"  yielding  a  very  large  gain  much  larger  than  d. 
Furthermore,  the  two  outcomes  are  equally  likely,  i.e.,  each  has  a 
50-50  chance  of  happening.  Finally,  a  further  moderating  action  is 
available  that  would  have  the  effect  of  making  this  project  less  of 
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an  all-or-nothing  proposition.  In  particular,  if  the  project  would 
have  broken  even,  this  action  would  add  a  gain  of  one  unit  of  money 
(a  meaningfully  large  denomination)  of  present  value.  On  the  other 
hand,  if  the  project  would  have  made  it  big,  the  action  would  decrease 
the  very  large  present  value  by  units  of  money.  Should  management 
approve  this  moderating  action?  This  obviously  will  depend  on  the 
size  of  x1,  the  price  that  must  be  paid  for  this  "insurance."  The 
question  that  needs  to  be  answered  by  management  is,  "what  is  the 
break-even  value  of  x^,  below  which  the  action  should  be  approved 
and  above  which  it  should  not?"  Given  this  break-even  value,  call  it 
x*,  it  is  then  simple  to  verify  that  the  desired  estimate  of  b^  is 


To  estimate  b^,  management  should  consider  a  very  similar 
hypothetical  situation.  The  one  difference  is  that,  instead  of 
perhaps  breaking  even,  the  first  possibility  for  the  outcome  of  the 
project  is  that  a  very  large  loss  (d  or  larger)  in  present  value 
will  be  incurred.  Nevertheless,  the  firm  is  already  committed  to  this 
project.  However,  a  moderating  action  again  can  be  undertaken  now  if 
desired.  Thus,  If  the  very  large  loss  would  have  occurred,  this 
action  would  decrease  the  loss  by  one  unit  of  money.  As  before,  if 
the  very  large  gain  would  have  been  attained,  the  action  would  decrease 
the  gain  by  some  amount,  which  will  be  denoted  here  by  The 

corresponding  question  again  needs  to  be  answered  by  management, 
namely,  "what  is  the  break-even  value  of  x^,  below  which  the  moder¬ 
ating  action  should  be  approved  and  above  which  it  should  not?" 


Given  the  break-even  value  x*,  the  desired  estimate  becomes 


b2  =  X^±  ,  X»/x*  . 


The  second  model  for  U(p)  differs  from  the  above  one  only  in 
the  behavior  of  the  utility  function  as  p  grows  very  large  in  the 
negative  direction.  Thus,  rather  than  assuming  that  U(p)  converges 
to  a  linear  asymptote  +  bgP  as  this  happens,  it  is  instead 
assumed  that  U(p)  decreases  exponentially  as  p  -*  -  00 ,  as  shown 
in  Figure  5.  Therefore,  the  algebraic  form  of  the  function  is 


U(p)  -  ax  +  b3P  -  a^e 


where  the  constants  for  the  last  term  are  implied  by  the  convention 
that  U(0)  =  0  and  U^^o)  =  1.  (As  before,  it  is  assumed  that 
>  0  and  0  <  b1  <  1. )  It  is  then  straightforward  to  find  the 
first  three  derivatives  as 


U(l)(p)  =  bx  +  (1  -  bx)e  '  1  , 

U<2>(p,  .  -a  I***2  -'(  ^  1 


Uu;(p)  =  a. 


Since  the  marginal  utility  (or  marginal  negative  utility  for  p 
decreasing)  U^(p)  increases  indefinitely  as  p  rather  than 

approaching  a  constant,  this  model  exhibits  an  even  greater  aversion 
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to  very  large  losses  than  the  basic  model.  Therefore,  it  will  be 
referred  to  hereafter  as  the  "high  risk  aversion"  model. 

Just  as  for  the  basic  model,  the  high  risk  aversion  model  can 
also  be  expressed  in  terms  of  the  three  parameters,  d,  b  ,  and  b^, 
although  their  precise  meaning  is  somewhat  different  here.  From  a 
mathematical  viewpoint,  d  now  is  a  somewhat  arbitrary  large  constant 
(although  its  interpretation  from  management's  viewpoint  is  similar 
to  before),  b  is  the  slope  of  the  asymptote  a^  +  b^p  which  U(p) 
converges  to  as  p  -*  +  °°,  just  as  before,  but  b^  now  is  defined 
as 

b2  =  U(l)(-d)  , 

the  marginal  utility  at  a  loss  in  present  value  of  d.  Therefore, 

b2  a  bl  +  (1  ’  Ve 

so  that 


where  in  denotes  the  logarithm  to  the  base  e.  Letting  k  denote 


y 


a,  thereby  becomes 


k  = 


a  =  (1  -  b  )kd  . 


■A 


Hence,  the  asymptote  a^  +  t^p  must  intersect  the  45°  line  through 
the  origin  at  (kd,kd),  as  shown  in  Figure  5- 

From  the  viewpoint  of  management,  the  interpretation  of  the  judg¬ 
ment  decisions  required  in  order  to  estimate  d,  b^,  and  is 

essentially  the  same  as  for  the  basic  model.  Thus,  d  may  be  inter¬ 
preted  to  management  just  as  before,  especially  if  it  is  not  clear  at 
the  outset  which  of  the  two  models  is  appropriate.  If  it  is  already 
apparent  that  the  high  risk  aversion  model  should  be  used,  then  there 
is  some  additional  flexibility  in  how  the  question  may  be  posed.  For 
example,  d  may  be  interpreted  simply  as  a  "very  large  loss"  above 
which  the  situation  would  be  deteriorating  disastrously.  Since  d 
is  only  a  convenient  benchmark  value  for  this  model,  the  important 
consideration  is  that  it  represent  a  magnitude  of  loss  that  is  mean¬ 
ingful  to  management  for  measuring  its  aversion  to  risk  (in  the  process 
of  estimating  bQ).  The  approach  to  estimating  bn  should  be  precisely 
the  same  as  for  the  basic  model.  The  approach  for  b^  also  should 
be  the  same  with  the  one  exception  that  the  unfavorable  possible  outcome 
from  the  hypothetical  project  must  now  be  a  loss  of  d,  rather  than 
any  very  large  loss  (d  or  larger)  as  before. 

Since  the  difference  between  the  two  models  is  in  the  behavior  of 
U''~'(p)  (which  can  also  be  interpreted  as  the  rate  of  increase  of 
negative  utility  when  p  is  decreasing)  as  p  -*  -  »,  the  choice 
between  them  would  be  based  on  an  analysis  of  this  behavior.  Thus, 
if  there  exists  a  real  disaster  level  where  further  losses  cannot  be 
absorbed  reasonably,  so  that  increasing  the  loss  above  this  level 
increases  the  damage  (negative  utility)  at  an  exponential  rate,  then 
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the  high  risk  aversion  model  should  be  used.  On  the  other  hand,  if 
very  large  possible  losses  could  be  absorbed  (albeit  reluctantly),  so 
that  increasing  such  a  loss  even  further  would  increase  the  detrimental 
effect  at  a  nearly  constant  rate,  then  the  basic  model  should  be  used. 
This  would  presumedly  be  the  case  if  the  total  commitment  for  the 
proposed  investments  would  be  quite  small  relative  to  the  corporate 
resources.  .  If  it  is  not  clear  which  case  applies,  then  management’s 
attitude  can  be  assessed  in  the  process  of  obtaining  the  judgment 
decision  leading  to  the  estimate  of  b^.  In  particular,  after  the 
hypothetical  situation  has  been  posed  in  terms  of  the  two  possible 
outcomes  (a  very  large  -loss  d  or  a  very  large  gain),  and  the  decision 
on  has  been  reached,  then  management  should  be  asked  if  their 

answer  would  change  much  if.  the  size  of  the  loss  for  the  first  possible 
outcome  were  increased  significantly.  An  answer  of  yes  would  suggest 
that  u(l}(p)  is  continuing  to  grow  substantially  as  p  grows  very 
large  in  the  negative  direction,  so  the  high  risk  aversion  model  should 
be  chosen.  An  answer  of  no  would  suggest  that  IJ^(p)  is  approaching 
a  constant  as  .  p  -*  -  °°,  so  the  basic  model  would  be  the  appropriate 
one . 

4.  Finding  the  Best  Combination  of  Investments 

Section  2  developed  a  model  for  considering  interrelationships 
between  investments  that  allows  one  to  identify  the  probability  dis¬ 
tribution  of  present  value  P(5)  for  given  feasible  solutions  5. 

In  particular,  expressions  were  obtained  for  the  mean  |i(s)  and 
2  '  -  . 

variance  a  (8)  of  this  distribution,  and  it  was  indicated  that  the 
functional  form  of  the  distribution  would  be  at  least  approximately 
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normal  under  certain  conditions.  Section  3  then  described  how  to 
construct  an  appropriate  utility  function  u(p)  (as  a  function  of  " i. 

the  realized  present  value  p)  which  formally  quantifies  management's 
judgment  on  risk  tradeoffs.  This  reduced  the  problem  under  consider¬ 
ation  to  that  of  determining  which  feasible  combination  of  investments 
maximizes  expected  utility,  i.e., 

max  E(U(P(5))) 

over  all  feasible  solutions  5.  What  remains  now  is  to  indicate  how 
E(u(P(8))}  can  be  calculated  for  a  given  8,  and  then  what  solution 
procedures  can  be  used  to  find  the  particular  feasible  8  which 
maximizes  this  quantity. 

Under  most  circumstances,  finding  the  exact  value  of  E(U(P(&))} 
is  not  a  straightforward  task,  since  this  requires  calculating  the 
expected  value  of  a  rather  complicated  function  u(*)  of  a  random 
variable  P(8)  having  perhaps  a  very  complex  probability  distribution. 
However,  one  exception  is  when  P(8)  has  a  normal  distribution  and 
u(p)  has  the  form  described  by  the  high  risk  aversion  model.  For 
this  case,  the  expected  utility  is  found  in  the  companion  volume 
[9>  P-  *0]  to  be 


E(u(P(8) ) )  =  ax  +  b1d(8)  -  a^e 

where  a  and  b^  are  defined  in  the  preceding  section.  In  most 
other  cases,  it  is  necessary  to  develop  an  approximation  of  E(U(P(5))} 
As  shown  in  [9>  PP-  35-31])  this  can  be  done  very  easily  by  using  a 


f  1-b, 


1-b,\2 

n(8)+  ~J  °  (8)/2 
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Taylor  series  expansion  and  ignoring  terms  above  the  second-order,— ^ 
which  yields 

U(2)(p(6» 

E{U(P(6))}  *U(n(5))  +  - g -  (&)  • 

This  approximation  is  a  very  convenient  one  since  it  only  requires 

knowing  the  mean  and  variance  of  P(&),  the  utility  function  U(p), 

(2) 

and  its  second  derivative  lr  '(p),  all  of  which  have  already  been 
identified  for  the  models  considered  here.  The  effect  of  using  the 
approximation  is  investigated  at  the  end  of  the  next  section. 

Two  solution  procedures  for  seeking  the  optimal  solution  (i.e., 
the  8  corresponding  to  the  best  combination  of  investments)  with 
respect  to  the  model  are  developed  in  the  companion  volume  [9>  pp.  49- 
67]*  Since  both  procedures  are  designed  to  be  executed  on  a  digital 
computer,  only  a  brief  introduction  to  their  general  nature  will  be 
repeated  here. 

One  of  these  methods  is  an  approximate  procedure,  i.e.,  it  will 
find  a  solution  8  which  should  be  a  very  good  one  but  is  not  guaranteed 
to  be  actually  optimal.  It  involves  solving  a  sequence  of  five  linear 
programming  problems,  where  each  one  is  constructed  by  introducing  a 
suitable  linear  approximation  to  the  objective  function  (e{u(P(S)))) 
based  on  the  solution  to  the  preceding  problem.  For  two  of  these 
problems,  it  is  necessary  to  obtain  an  approximate  integer  solution, 
which  can  be  done  either  by  rounding  the  linear  programming  solution 
or  by  using  suitable  integer  programming  methods.  (These  methods 

—  See  [9>  P-  36]  for  this  remainder  term  expressed  in  terms  of  the 

t 

third  derivative  Uw'(p). 
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include  the  highly  efficient  heuristic  procedures  developed  elsewhere 

integer  solutions  (usually  the  second  one)  is  the  desired  approximate 
solution. 

Since  computer  codes  for  solving  a  linear  programming  problem  are 
widely  available  and  very  efficient,  the  approximate  procedure  provides 
a  convenient  approach  that  can  even  be  used  on  extremely  large  problems. 
Furthermore,  as  Weingartner  [17]  has  found,  linear  programming  is  a 
powerful  tool  for  the  further  analysis  of  the  problem.  For  example, 
the  dual  evaluators  provide  vital  information  for  studying  whether  ib 
would  be  worthwhile  to  increase  the  initial  allocation  of  scarce 
resources  (such  as  budgeted  funds)  available  for  the  proposed  invest¬ 
ments.  Another  crucial  advantage  is  the  ease  with  which  an  extensive 
sensitivity  analysis  of  the  input  parameters  of  the  model  may  be 
conducted. 

The  second  solution  procedure  is  an  exact  one,  i.e.,  the  solution 

5  it  obtains  is  definitely  optimal  with  respect  to  the  model  and 

objective  function  used.  It  is  based  on  the  branch-and-bound 
17/ 

technique, —  inasmuch  as  the  procedure  consists  of  applying  a  special 

"l8/ 

new  branch-and-bound  algorithm  for  integer  nonlinear  programming — 
to  first  a  subproblem  and  then  to  the  original  problem.  The  objective 
of  the  subproblem  is  to  find  the  feasible  solution  6  which  maximizes 

___  ■ 

—  Gee  Hlllier  and  Lieberman  [11,  pp.  565-570]  for  an  elementary 
description  of  this  technique. 

—  This  algorithm  is  presented  in  the  companion  volume  [9,  pp-  53-60]. 


b^^he^uthor^[6j^or,  if  the  problem  is  not  too  large,  the  algorithm 
given  by  Petersen  [15]  and  Geoffrion  [  3  ]  • )  "T^nTTeETer" 
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expected  present  value  E{P(8)}.  This  solution  is  used  to  structure 
the  original  problem  in  such  a  way  that  the  branch-and-bound  algorithm 
can  then  be  applied  directly  to  it. 

Although  the  exact  procedure  is  not  computationally  feasible  for 
the  very  large  problems  that  can  be  handled  by  the  approximate  pro¬ 
cedure,  it  is  quite  efficient  for  problems  of  moderate  size  (approxi¬ 
mately  ten  investment  proposals).  The  branch-and-bound  approach  is 
not  as  well-suited  as  linear  programming  for  sensitivity  analysis, 
etc.,  so  it  may  sometimes  be  advantageous  to  apply  both  procedures. 
Furthermore,  the  two  integer  solutions  obtained  in  the  course  of  the 
approximate  procedure  can  be  used  to  expedite  the  respective  appli¬ 
cations  cf  the  branch-and-bound  algorithm  in  the  exact  method. 

Computational  experimentation  measuring  the  efficiency  of  the 
two  solution  procedures,  and  the  effectiveness  of  the  approximate 
procedure  in  obtaining  good  solutions,  is  reported  in  the  next  section. 

5-  Computational  Experience 

In  order  to  evaluate  the  two  solution  procedures  described  in  the 
preceding  section,  ALGOL  codes  were  developed  for  the  IBM-36o/67«— 1 ' 

For  the  two  phases  of  the  approximate  procedure  where  an  approximate 
integer  solution  is  required,  one  of  the  heuristic  procedures  for 
integer  programming  developed  by  the  author  [6]  was  used. — The 


■ 

—  The  author  wishes  to  express  his  great  appreciation  to  Ilan  Adler 
and  Mrs.  Patricia  Peterson  for  developing  the  codes,  and  to  Francisco 
Pereira  for  running  them  to  obtain  the  computational  results  given  in 
this  section. 


20/ 


The  procedure  used  is  labeled  1-2A-1  in  [6]. 
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main  objective  in  designing  the  codes  was  reliability  rather  than 
efficiency  or  operational  convenience. 

The  next  step  was  to  generate  all  of  the  parameters  described  in 
Section  2  for  ten  basic  test  problems,  each  involving  ten  investment 
proposals  (m  =  lo) .  This  was  done  on  the  computer  by  randomly 
generating  some  integers  from  each  of  the  intervals  indicated  in 

2i  / 

Table  I. —  In  the  case  of  the  p.,  and  p.,  (k  >  j),  some  of 

jK-  jK 

them  were  randomly  selected  to  have  a  value  of  zero  according  to  a 

prespecified  probability  given  in  the  second  column  of  Table  I,  and 

only  the  rest  of  them  were  assigned  randomly  generated  integers.  In 

2 

the  case  of  the  a . ,  another  parameter  s,  was  generated  instead, 

1  J 

2 

and  then  used  to  calculate  a  in  different  ways  to  obtain  a  variety 

J 

of  risk  profiles,  as  will  be  described  shortly.  Each  of  the  test 

m 

problems  includes  five  constraints  of  the  form,  J)  a  6  <  b 

fcTi  K  J 

(which  might,  for  example,  represent  budget  constraints  in  the  first 

five  time  periods),  and  these  a.,  and  b,  also  were  randomly 

dK  J 

generated  from  the  intervals  indicated  in  Table  I. 


—  A  square  bracket  at  the  end  of  an  interval  indicates  that  the 
end-point  is  included  among  the  integers  eligible  for  random  selection, 
whereas  a  round  parenthesis  indicates  that  it  is  not  included. 
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Table  I 


Description  of  the  Randomly  Generated 
Parameters  for  the  Basic  Test  Problems 


Type  of 

Prob,  of  Presetting 

Interval  Used  for 
Random  Generation 

^  JL  Kt  0  U  U 

0 

[-10,90) 

^k 

1/2 

(-25,25) 

Sj 

0 

[0,100) 

10pjk 

1/2 

(-10,10) 

alk 

0 

[0,100) 

n2k 

0 

[-10,90) 

a3k 

0 

[-20,80) 

aUk 

0 

[-30,70) 

a5k 

0 

[-40, 60) 

0 

[125,375) 

In  order  to  fully  define  the  problems ,  it  was  next  necessary  to 
2 

calculate  the  a  in  terms  of  the  s  and  to  set  the  values  of  the 

J  J 

three  utility  function  parameters  described  in  Section  5.  This  was 
done  in  three  different  ways,  as  shown  in  Table  II,  in  order  to 
obtain  sets  of  problems  involving  investment  proposals  with  very 
little  risk,  a  moderate  degree  of  risk,  and  a  very  high  level  of 
risk  associated  with  them,  respectively.  The  indicated  value  of 
by  (which  is  i  ^  e  rounded  to  three  significant  digits)  was  so 
chosen  in  order  to  have  k  =  1  in  the  high  risk  aversion  model  for 
U(p).  This  yields  the  same  linear  asymptote  when  p  -*  «  as  for  the 
basic  model  for  U(p) ,  thereby  providing  a  better  comparison  of 
other  features  of  the  two  models. 


Table  II 

Value  of  the  Fixed  Parameters  for  the 
Basic  Test  Problems  % 


Parameter 

Low  Risk 
Problems 

Moderate  Risk 
Problems 

High  Risk 
Problems 

2 

a. 

3 

a2=s. 

J  J 

<?  =  ,?. 

d 

4o 

ho 

200 

V 

0.5 

0.5 

0.5 

b2 

1.86 

1.86 

1.86 

The  resulting  30  problems  were  then  used  to  test  the  effectiveness 
of  the  approximate  procedure  by  comparing  its  solution  with  that 
obtained  from  the  exact  procedure,  as  well  as  with  the  exact  solution 
for  maximizing  p(5)  if  risk  were  ignored  (which  is  obtained  from 
the  first  phase  of  the  exact  procedure).  This  was  done  first  with 
the  basic  model  for  U(p),  and  then  with  the  high  risk  aversion  model 
for  U(p)  when  the  distribution  of  P(8)  is  assumed  to  be  normal. 

The  results  are  shown  in  Tables  III,  IV,  and  V  for  the  three  types  of 
problems.  For  each  protlem  and  approach,  the  table  gives  (l)  "E," 
the  value  of  E(U(P(5))}  for  the  solution  obtained,  (2)  "p,"  the 
value  of  p(5)  for  this  solution,  (3)  "d,"  the  value  of  0(5)  for 
this  solution,  and  (4)  "Time,"  the  time  in  seconds  required  to  obtain 
the  solution  on  an  IBM  360/67  computer.  For  each  model,  the  three 
columns  in  the  tables  refer  respectively  to  (l)  the  feasible  solution 
which  maximizes  p(o)  without  considering  risk,  (2)  the  solution 
obtained  by  the  approximate  procedure,  and  (3)  the  optimal  solution 


as  obtained  by  the  exact  procedure. 


Table  III 


Performance  of  the  Approximate  Procedure 
with  Low  Risk  Test  Problems 


ft 

Basic  ’ 

Basic  Model 

High  Risk  Aversion  Model 

HU 

Approx. 

Proc. 

Exact 

Proc. 

max  n(8) 

Approx. 

Proc. 

Exact 

Proc. 

#1 

E 

h 

a 

Time 

213.0 

1*25.0 

10.6 

0.73 

212.5 

424,0 

8.8 

O.78 

213.1 

425.0 

10.6 

3*02 

232.5 

425.0 

10.6 

0-73 

232.0 

424.0 

8.8 

0.70 

232.5 

425.0 

10.6 

3.04 

#2 

E 

P 

a 

Time 

26k. 2 
528.0 
22.9 
0.59 

237.9 

475.0 

15.7 

0.64 

264.2 

528,0 

22.9 

1.20 

284.0 

528.0 

22.9 

0.59 

257.5 

475.0 

15.7 

0.6l 

284.0 

528.0 

22.9 

1.47 

#3 

E 

P 

a 

Time 

142.5 

284.0 

6.6 

1.38 

87.7 

176.0 

20.5 

0.86 

142.5 

284.0 

6.6 

8.01 

162.0 

284.0 

6.6 

1.38 

107-7 

176.O 

20.5 

0.85 

162.0 

284.0 

6.6 

7.74 

#4 

E 

b 

0 

Time 

152.9 

305.0 

12.9 

l.lo 

152.9 

305.0 

12.9 

0.64 

152.9 

305.0 

12.9 

3.08 

172.5 

305.0 

12  9 
1.10 

172.5 

305-0 

12.9 

0.66 

172.5 

305.0 

12.9 

3.69 

#5 

E 

P 

c 

Time 

3:9.5 

259.0 

19.8 

0.93 

128.9 

258.0 

21.9 

0.67 

129.5 

259.0 

19,8 

2.70 

149-5 

259-0 

19.8 

0.93 

149.0 

258.0 

21.9 

0.75 

149-5 

259.0 

19.8 

2.67 

#6 

E 

b 

0 

Time 

1B1- 

Wttm 

112.4 

225.0 

21,0 

0.75 

117.1 

234.0 

16.9 

7-65 

136.9 

234.0 

17.8 

1.63 

132.4 

225,0 

21.0 

0-77 

137.0 

234.0 

16.9 

7.64 

#7 

E 

r 

o 

Time 

53.1 

lo6.o 

12.0 

0.78 

53.1 

lu6.0 

12,0 

0.75 

53,1 

106,0 

12.0 

1-95 

71.5 

106.0 

12.0 

O.78 

71.5 
106.0 
12.0 
0.6  7 

71.5 

106.0 

12.0 

1.99 

:-8 

E 

P 

a 

T  ime 

167.6 

335.0 

20.9 

0.85 

167,6 
;35  ■  0 

20,9 

0.68 

I67.6 

335.0 

20.9 

4.30 

187.5 

335.0 

20.9 

0.85 

187.5 
335.  c 
20.9. 
0.66 

167.5 
335-0 
20 .9 
4.17 

I 


j 


E 

190.8 

190.8 

190.8 

210.5 

210.5 

210.5 

#9 

H 

381.0 

381.0 

381.0 

381.0 

381.0 

381.0 

Cf 

16.4 

16.4 

16.4 

16.4 

16.4 

16.4 

Time 

2.83 

0.71 

8.10 

2.83 

0.74 

8.03 

E 

128.5 

121.7 

128.5 

148.0 

mmm 

148.0 

#10 

H 

256.0 

243.0 

256.0 

256.0 

256.0 

a 

8.7 

15.9 

8.7 

8.7 

8.7 

Time 

0.78 

0.71 

4.24 

0.78 

4.21 

E 

155-9 

146.6 

156.0 

175.5 

166.2 

175.5 

Average 

H 

311.3 

292.8 

311.3 

311.3 

292.8 

311-3 

a 

14.9 

16.6 

14.8 

14.9 

16.6 

14.8 

Time 

1.16 

0.72 

3.42 

0.71 

4.46 

Table  IV 


Performance  of  the  Approximate  Procedure 
v-ith  Moderate  Risk  Test  Problems 


Basic 

Basic  Model 

High  Risk  Aversion  Model 

Probl.  . 

Approx. 

Exact 

Approx. 

Exact 

Used 

max  |i(6) 

Proc. 

Proc. 

max  p(6) 

Proc. 

Proc. 

E 

207-9 

208.0 

208.0 

232.5 

232.5 

232.5 

#1 

h 

425.0 

424.0 

424.0 

425.0 

425.0 

425.0 

a 

80.2 

75.4 

75.4 

80.2 

80.2 

80.2 

Time 

0.73 

0.73 

2.71 

0.73 

0.76 

3.6l 

E 

242.9 

228.5 

242.9 

282.6 

257.5 

282.6 

#2 

528.0 

475.0 

528.0 

528.0 

475.0 

528.0 

a 

183.5 

115.9 

183-5 

183.5 

115.9 

183.5 

Time 

0.59 

0.65 

2.46 

0.59 

0.6l 

2.19 

E 

139.5 

111.9 

139-5 

162.0 

-109-6 

162.0 

#3 

H 

284.0 

236.0 

284.0 

284.0 

176.0 

284.0 

a 

50.5 

67.8 

50.5 

50.5 

147.3 

50.5 

Time 

1.58 

O.78 

.  7.4o 

1.38 

0.85 

7-99 

E 

142.0 

144.3 

144.3 

172.3 

172.3 

172.3 

#4 

305.0 

289.0 

289.0 

305.0 

305.0 

305.0 

cr 

99.6 

24.8 

24.8 

99.6 

99  -6 

99-6 

Time 

1.10 

0.66 

3.27 

1.10 

0.66 

3.15 

E 

95.6 

97-0 

110.9 

34.7 

-724.6 

134.3 

#5 

h 

259-0 

216.0 

22Q.0 

259.0 

258.0 

229.0 

a 

162.2 

85.7 

53.1 

162.2 

180.9 

53.1 

Time 

0.93 

0.68 

2.88 

0.93 

0.67 

3-75 

E 

93-0 

93-0 

101.5 

125.5 

-522.4 

133.4 

#6 

h 

254.0 

234.0 

229.0 

234.0 

225.0 

229.0 

a 

130.2 

130.2 

95-5 

130.2 

170.8 

95.5 

Time 

1.63 

0,76 

6.25 

1.63 

0-79 

8.05 

E 

20.9 

30.1 

30.1 

38.5 

57-5 

57-5 

#7 

h 

106.0 

75-0 

75.0 

106.0 

98.0 

98.0 

0 

101.1 

42.0 

42.0 

101.1 

78.0 

78.0 

Time 

O.78 

0.60 

1.82 

O.78 

O.69 

2.13 

E 

140.7 

l4o.7 

140.7 

166.2 

166.2 

166.2 

#8 

|i 

335.0 

335.0 

335.0 

335.0 

335.0 

335.0 

0 

164.3 

164.3 

164.3 

164.3 

164.3 

164.3 

Time 

0.85 

O.69 

4.72 

0.85 

0.71 

5.01 

E 

178.6 

178.6 

178.6 

210.4 

210.4 

210.4 

#9 

381.0 

381.0 

381.0 

381.0 

381.0 

381.0 

a 

118.4 

118.4 

118.4 

118.4 

118.4 

118.4 

Time 

2.83 

0.70 

7.35 

2.83 

0.72 

8.03 

E 

124.6 

112.9 

124.6 

147.9 

134.8 

147-9 

#10 

H 

256.0 

241.0 

256.0 

256.0 

243.0 

256.0 

a 

54.9 

76.0 

54.9 

54.9 

126.2 

54.9 

Time 

0.78 

0.69 

4.25 

0.78 

0.72 

4.36 

E 

138.6 

134.5 

142.1 

157.2 

-12.5 

169.9 

Average 

H 

277.8 

290.6 

221.1 

277-8 

292.1 

307.0 

a 

114.5 

90.0 

86.2 

114.5 

128.2 

97-8 

Time 

1.16 

0.69 

4.31 

1.16 

0.72 

4-84 

Table  V 


Performance  of  the  Approximate  Procedure 
with  High  Risk  Test  Problems 


Basic 

Problem 

Used 

Basic  Model 

High  Risk  Aversion  Model 

max  n(&) 

Approx. 

Proc. 

Exact 

Proc. 

max  u(B) 

Approx. 

Proc. 

Exact 

Proc. 

#1 

E 

a 

Time 

-101 

425 

630 

0.73 

40.0 

95 

53-8 

O.76 

56.9 

28l 

268 

3.65 

-1384 

425 

630 

0.73 

-1082 

424 

617 

0.73 

180.3 

281 

268 

6.60 

#2 

E 

H 

a 

Time 

-1050 

528 

1431 

0.59 

-7558 

370 

1017 

0.64 

0 

0 

0 

4.91 

-9.8X1011 

528 

1431 

0.59 

-734.6 

314 

556 

O.60 

64.6 

91. 

139 

24.98 

#3 

E 

a 

Time 

-112 

284 

465 

1.38 

-212 

236 

482 

0.82 

77-9 

219 

149 

6.92 

-118 

284 

465 

1.38 

-343 

236 

482 

O.78 

166.6 

276 

289 

6.00 

#4 

E 

a 

Time 

-541 

305 

792 

1.10 

n 

Warn 

mm 

-55,571 

305 

792 

1.10 

-35,230 

231 

750 

0.71 

164.5 

261 

267 

3.43 

#5 

E 

a 

Time 

-2468 

259 

14 11 

0.93 

48.2 

118 

68.9 

0.60 

48.2 

118 

68.9 

5.74 

-6.2X1010 

259 

14 11 

0.93 

100.2 

118 

68.9 

0.6l 

100.2 

118 

68.9 

5.96 

#6 

E 

H 

a 

Time 

-1207 

234 

959 

1.63 

27.3 

51 

15.6 

0.63 

27.3 

51 

15.6 

4.53 

-3,042,622 

234 

959 

1.63 

19 

lEHlB 

52.6 

81 

152 

5.46 

#7 

E 

H 

a 

Time 

-2261 

106 

856 

0.78 

-47.2 

43 

96.2 

0.59 

0 

0 

0 

1.61 

-554,102 

106 

856 

0.78 

-36.0 

75 

272 

0.66 

30.9 

43 

96.2 

1.70 

#8 

HI 

-1622 

335 

1332 

0.85 

48.6 

142 

103 

0.67 

48.6 

142 

103 

2.70 

-8.1X1010 

335 

1332 

0.85 

114.9 

142 

103 

0.65 

114.9 

142 

103 

3.32 

48 


E 

-514 

-473 

51.2 

-312,930 

-239 

160.6 

#9 

H 

38l 

122 

146 

381 

250 

269 

a 

892 

442 

103 

892 

472 

289 

Time 

2.83 

0.64 

8.91 

2.83 

0.67 

9.45 

E 

-60.5 

95.7 

95-7 

57-9 

57.9 

154.7 

#10 

256 

188 

188 

256 

256 

188 

a 

381 

24.0 

24.0 

381 

381 

24.0 

Time 

0.78 

0.64 

3.53 

0.78 

0.70 

3.35 

E 

-995.6 

-176.9 

45.7 

-1.123X1011 

.37344.4 

119.0 

Average 

311.3 

167.0 

140.6 

311.3 

209.7 

175.0 

a 

914.9 

309.5 

99.9 

914.9 

371.8 

159-6 

Time 

1.16 

0.60 

4.26 

1.16 

0.68 

7.51 

49 


An  examination  of  these  three  tables  reveals  that  the  approximate 
procedure  obtained  the  optimal  solution  in  slightly  less  than  half  of 
the  problems.  For  the  low  risk  test  problems,  optimality  was  attained 
in  four  of  the  ten  cases  with  both  models  for  U(p)>  and  the  objective 
function  was  within  1C$  of  the  optimal  value  in  five  of  the  remaining 
six  problems.  Hie  optimal  solution  was  obtained  for  five  of  the  ten 
moderate  risk  problems.  For  the  high  risk  test  problems,  the  optimal 
solution  for  the  basic  model  was  obtained  four  times,  and  the  optimal 
solution  for  the  high  risk  aversion  model  was  obtained  in  two  cases, 
but  the  solution  was  far  from  optimal  for  most  of  the  other  problems 
of  this  kind.  However,  the  approximate  procedure  was  successful  in 
surpassing  the  max  u(5)  solution  (the  best  solution  if  risk  were 
ignored)  in  eight  of  the  ten  high  risk  test  problems  (with  a  tie  in 
one  of  the  remaining  two  cases). 

The  results  described  are  all  for  problems  having  essentially  the 
same  set  of  values  for  the  parameters  of  U(p).  Therefore,  a  number 
of  different  changes  were  made  in  these  parameter  values  and  incorpor- 
ated  into  the  moderate  risk  version  (i.e.,  o_.  =  ^ )  of  test  problem 
5  and  the  high  risk  version  (i.e.,  o ^  =  s^)  of  test  problem  6.  The 
objective  was  to  gain  some  insight  into  both  the  effectiveness  of  the 
approximate  procedure  and  the  behavior  of  the  solutions  as  the  utility 
function  changes.  The  results  are  shown  in  Tables  VI  and  VII,  where 
the  first  three  columns  indicate  the  parameter  values  used  for  each 
run.  In  order  to  provide  a  convenient  standard  of  comparison,  the 
first  set  of  rows  repeats  the  results  with  the  original  set  of  parameter 
values  from  the  preceding  tables. 


Table  VII 


Effect  of  Changing  the  Utility  Function 
on  High  Risk  Test  Problem  6 


l 

|  Basic  Model  f 

|  High  Risk  Aversion  Model 

d  bJ 

ba 

■ 

Approx. 

Proc. 

max  n(8) 

Approx. 

Proc. 

Exact 

Proc. 

200  0.5 

1.86 

n 

-1207 

234 

959 

1.63 

27.3 

51 

15.6 

0.63 

27.3 

51 

15.6 

4.53 

-3,042,622 

234 

959 

1.63 

47.9 

51 

15.6 

0.64 

52.6 

81 

152 

5.46 

! 

100  0.5 

1.86 

E 

a 

Time 

-1213 

234 

959 

1.63 

26.0 

51 

11.2 

0.66 

26.0 

51 

11.2 

4.54 

-4.4x10*° 

234 

959 

1.63 

45.3 

51 

11.2 

O.67 

45.3 

51 

11.2 

5-87 

400  0.5 

1.86 

E 

H 

a 

Time 

-1197 

234 

959 

I.63 

-880 

105 

545 

0.65 

29.6 

51 

11.2 

4.55 

-1655 

234 

959 

1.63 

65.5 

81 

149 

0.66 

84.3 

151 

326 

5.04 

200  0.1 

1.86 

E 

H 

a 

Time 

-1709 

234 

959 

1.63 

30.6 

51 

11.2 

0.65 

30.6 

51 

11.2 

4.43 

-21,108 

234 

959 

1.63 

47.1 

51 

11.2 
0.66 

47.1 
51 

11.2 
5-04 

o\ 

6 

8 

CJ 

1.86 

E 

h 

a 

Time 

-1103 

234 

959 

1.63 

-1584 

105 

545 

0.66 

n 

■ 

-2.  cxio^5 
234 

959 

1.63 

-5-ixio^ 

225 

533 

0.67 

>  180 

200  0.5 

1.11 

E 

a 

Time 

-84l 

234 

959 

1.63 

48.3 

51 

11.2 

0.72 

48.3 

51 

11.2 

4.79 

-10.6 

234 

959 

1.63 

74.3 

81 

149 

0.66 

122.0 

151 

326 

4.30 

200  0.5 

10 

E 

H 

a 

Time 

-9551 

234 

959 

1.63 

13-7 

51 

11.2 

0.64 

13.7 

51 

11.2 

4.4l 

-2.CK10  J 
234 

959 

1.63 

-38.3 

81 

149 

O.65 

43.2 
51 

11.2 
785 

The  next  objective  was  to  test  the  efficiency  of  the  exact  pro¬ 
cedure  with  problems  of  different  sizes.  It  was  seen  in  the  preceding 
tables  that  it  is  quite  efficient  with  problems  having  five  constraints 

and  ten  variables,  requiring  between  one  and  ten  seconds  in  all  but 

22/ 

one  of  the  84  cases. — '  However,  the  nature  of  this  branch-and- bound 
procedure  is  that. its  running  time  should  tend  to  grow  quite  rapidly 
with  the  size  of  the  ^problem  (especially  as  m,  the  number  of  variables 
increases).  In  order  to  estimate  this  rate  of  growth,  and  the  resulting 
maximum  number  of  variables  (investment  proposals)  for  which  the  pro¬ 
cedure  should  be  computationally  feasible,  some  of  the  basic  test 
problems  were  modified  and/or  combined  to  construct  problems  of  other 
sizes.  In  each  case,  the  values  of  the  parameters  shown  in  Table  II 
were  retained  with  the  one  exception  that  d  was  multiplied  by  -jL 

All  of  the  original  values  of  the  p  ,  p  ,  s  ,  and  p  also  were 

0  JK  J  JK 

retained.  When  new  values  of  p^k  and  had  to  be  introduced 

because  of  combining  problems,  these  new  values  were  automatically 
set  equal  to  zero.  All  variables  kept  their  same  a  coefficients 

JK 

in  the  constraints.  In  those  cases  were  the  number  of  constraints 
was  reduced  from  five  to  one,  this  was  done  by  eliminating  all  but  the 
constraint  involving  the  a^  and  b^  parameters.  When  a  basic  test 
problem  was  used  to  contribute  only  half  of  its  variables  to  the  new 
problem,  the  first  five  variables  were  chosen,  and  the  values  of  the 
b^  were  divided  by  two.  When  combining  two  basic  test  problems  or 

___ 

— -  When  the  high  risk  aversion  model  was  used  with  the  b^  =  0.9 

problem  in  Table  VII,  the  procedure  did  not  terminate  in  three  minutes 
of  computer  time.  The  reason  is  not  evident. 
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portions  thereof,  their  current  values  of  b^  were  added  to  obtain 
the  ^  for  the  new  problem. 

Table  VIII  shows  the  results  of  the  runs  on  an  IIM-360/67  system 
with  these  new  test  problems.  The  first  column  indicates  which  basic 
test  problea(s)  was  used  to  construct  the  problem  for  that  run,  tod 
the  next  two  columns  give  the  resulting  number  of  constraints  and 
variables.  Th2  remaining  three  pairs  of  columns  identify  the  version 
of  the  problem  (corresponding  to  the  three  versions  in  Table  II)  and 
the  type  of  utility  function  model  (basic  model  or  high  risk  aversion 
model)  that  were  used.  To  provide  a  frame  of  reference,  the  average 
time  for  all  of  the  basic  test  problems,  as  well  as  the  individual 
times  for  the  ones  used  to  construct  new  problems  here,  are  repeated 
from  Tables  III,  IV,  and  V. 

Recall  from  the  preceding  section  that  an  exact  expression  for 
E(U(P(b))}  is  now  known  only  when  the  high  risk  aversion  model  is 
being  used  and  the  distribution  of  P(8)  is  normal.  Otherwise,  only 
an  approximate  expression,  obtained  from  using  a  Taylor  series  expansion, 
is  available.  The  exact  procedure  finds  the  feasible  solution  that 
maximizes  whatever  objective  function  is  being  used.  Therefore,  when 
the  approximate  expression  is  taken  to  be  the  objective  function,  there 
is  no  guarantee  that  the  resulting  solution  will  be  the  exact  optimal 
solution  for  the  real  problem.  Thus,  it  is  important  to  verify  that 
this  approximation  really  is  an  adequate  one,  especially  in  terras  of 
the  resulting  solution  still  being  at  least  nearly  optimal. 

It  is  straightforward  to  study  this  question  analytically  for 
the  case  of  the  high  risk  aversion  model.  Define 
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Time  (in  sec.)  Required  by  the  Exact  Procedure 
for  Test  Problems  of  Different  Sizes 
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Thus,  the  approximate  expression  is  equivalent  here  to  expanding  the 
c 

e  part  of  the  exact  expression  and  then  ignoring  terms  above  the 
first-order.  The  effect  is  to  give  somewhat  less  weight  to  the  penalty 
for  risk  (as  measured  by  cj(6))  than  is  warranted.  The  approximation 
is  quite  good  if  c  is  somewhat  less  than  one,  but  is  quite  poor  if 

c  is  somewhat  larger  than  one.  Thus,  the  critical  factor  is  the 

a  a. 

relative  size  of  ^Tb~  and  the  Recall  that-  =  kd,  and 

that  (kd,kd)  is  the  point  at  which  the  asymptote  (a1  +  b^p)  of  the 
utility  function  intersects  the  ^5°  line  through  the  origin.  Thus, 

al 

j  ^  ■  can  be  interpreted  loosely  as  a  point  of  diminishing  returns  for 
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gains  in  present  value.  If  this  point  is  at  least  comparable  in  size 
to  0(6)  for  the  interesting  values  of  6,  then  the  approximation 
should  be  quite  adequate. 

For  the  basic  test  problems  and  high  risk  aversion  model,  the  value 
of  kd  is  4o  for  the  low  risk  and  moderate  risk  versions  and  200  for 
the  high  risk  version  of  these  problems.  The  corresponding  values  of 
the  c(b)  for  interesting  S  are  such  that  the  value  of  c  would  be 
a  rather  small  fraction  of  one  for  the  low  risk  problems,  but  would 
tend  to  be  approximately  one  or  even  considerably  larger  for  most  of 
the  other  problems.  Therefore,  if  the  approximate  expression  for 
E(U(P(8)))  were  to  be  used  instead  of  the  exact  one,  it  appears  that 
the  exact  procedure  should  still  give  very  good  solutions  for  the 
low  risk  problems  but  might  not  perform  well  for  the  other  problems. 

To  test  this  experimentally,  all  30  problems  were  rerun  with  the  exact 
procedure  using  the  approximate  expression.  The  results  are  shown  in 
Table  IX.  The  first  column  identifies  the  problem,  where  LR  is  an 
abbreviation  for  Low  Risk,  MR  for  Moderate  Risk,  and  HR  for  High 
Risk.  The  second  column  indicates  the  time  in  seconds  required  by  the 
exact  procedure  run  in. this  way,  and  the  next  two  columns  give  the 
values  of  |i(8)  and  0(8)  for  the  resulting  solution  which  maximizes 
the  approximate  objective  function.  The  fifth  column  recalculates 
the  value  of  the  objective  function  for  this  solution  by  using  the 
exact  expression.  The  last  column  then  gives  the  ratio  (expressed 
as  a  percentage)  of  this  value  to  the  corresponding  value  for  the 
real  optimal  solution  as  given  in  Tables  II,  III,  or  IV.  An  examination 
of  this  column  also  reveals  that  the  solution  obtained  actually  was 
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Table  IX 


optimal  for  all  ten  low  risk  problems,  and  for  seven  out  of  ten  of 
both  the  moderate  risk  and  high  risk  problems.  Studying  the  effect 
of  using  the  approximate  expression  with  the  basic  model  for  U(p) 
is  more  difficult.  However,  an  analysis  of  the  remainder  term  that 
has  been  deleted  from  this  expression  suggests  that  it  probably  would 
tend  to  be  positive,  and  should  tend  to  increase  as  cr(B)  increases, 
so  that  the  expression  underestimates  E{u(P(B))}  most  for  large  values 
of  cr(8).  Therefore,  it  appears  that,  in  contrast  to  the  conclusion 
drawn  for  the  high  risk  aversion  model,  the  effect  of  the  approximation 
with  the  basic  model  is  to  give  somewhat  more  weight  to  the  penalty 
for  risk  than  is  warranted.  This  conjecture  is  strongly  supported  by 
the  computational  results  given  in  this  section.  In  those  cases  where 
the  basic  model  (using  the  approximate  expression)  and  the  corresponding 
high  risk  aversion  model  (using  an  exact  expression)  obtained  different 
solutions  on  a  given  problem,  it  was  always  the  basic  model  that  yielded 
the  solution  with  the  smaller  value  of  u(B),  despite  the  less  con¬ 
servative  utility  function  for  this  model. 

6.  Conclusions 

The  companion  volume  l 9]  comprehensively  studies  the  problem  of 
evaluating  risky  interrelated  investments.  One  basic  approach  developed 
there  has  been  selected  for  further  study  and  elaboration  here.  This 
approach  involves  considering  the  probability  disti xbution  of  present 
value  for  alternative  combinations  of  investments,  and  then  using 
powerful  solution  procedures  to  seek  the  feasible  combination  that 
maximizes  the  expected  value  of  a  utility  function  of  present  value. 

It  takes  into  account  various  kinds  of  interrelationships  between 
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the  proposed  investments  in  some  relatively  convenient  ways. 

The  utility  function  approach  to  considering  risk  is  a  relatively 
sophisticated  technique  that  usually  has  been  shunned  in  practice  in 
favor  of  subjective  evaluation  of  risk  profiles  by  management.  However, 
when  a  substantial  number  of  proposed  projects  are  being  evaluated 
simultaneously  (due  to  such  interrelationships  as  competing  for  the 
same  capital  resources),  so  that  the  number  of  possible  investment 
combinations  is  huge,  it  is  then  necessary  to  use  some  other  device 
to  identify  the  few  best  alternatives  for  careful  evaluation  by  manage¬ 
ment.  The  utility  function  approach  is  a  powerful  and  flexible  way  of 
doing  this.  Furthermore,  a  practical  way  of  constructing  a  reasonable 
utility  function  by  merely  obtaining  three  meaningful  judgment  decisions 
from  management  has  been  developed  in  Section  Computational  experience 
(see  Tables  VI  and  VII)  also  suggests  that  the  optimal  solution  tends 
to  be  not  too  sensitive  to  the  choice  for  each  of  these  decisions. 
Therefore,  the  use  of  a  utility  function  in  this  way  appears  to  be 
a  much  more  attractive  and  practical  approach  to  the  problem  than 
might  have  been  anticipated. 

A  useful  area  for  future  research  would  be  a  behavioral  study  of 
the  process  by  which  management  would  make  the  three  judgment  decisions 
described  in  Section  5  and  of  hew  the  resulting  answers  would  relate 
to  management's  true  preferences  on  risk  tradeoffs.  This  might  lead 
to  a  refined  procedure  for  obtaining  the  required  information  from 
management.  » 

It  is  important  to  recognize  that,  when  trying  to  identify  the 
best  combination  of  investments,  it  is  the  risk  profile  for  the 
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combination  (primarily  p(0)  and  a (6) ).  that  is  crucial  rather  than 
the  risk  profiles  for  the  individual  investments.  Urns,  even  though 


a  number  of  proposals  appear  individually  to  be  very  risky  (i.e., 

is  large  relative  to  p.) ,  there  may  still  be  an  acceptably  small 

degree  of  risk  associated  with  approving  all  of  these  investments 

jointly.  This  diversification  effect  occurs  because  o(8)  tends  to 

23/ 

grow  less  rapidly  than  p(8)  as  more  investments  are  approved,— 
and  this  effect  is  even  more  pronounced  if  some  of  the  P.  are  nega- 

J 

tively  correlated.  Thus,  it  is  dangerous  to  rely  too  much  on  sub¬ 
jective  or  compartmentalized  decision-making  here. 

Computational  experience  suggests  that  the  approximate  solution 
is  fairly  effective  in  obtaining  good  solutions  (and  slightly  more  so 
when  using  the  basic  model  for  U(p)  rather  than  the  high  risk  aversion 
model),  since  it  actually  yielded  the  optimal  solution  for  almost 
half  of  the  test  problems.  However,  it  does  tend  to  give  somewhat 
less  weight  to  the  element  of  risk  than  is  warranted.  As  a  result, 
it  may  yield  a  poor  solution  when  considering  a  very  risky  set  of 
proposed  investments.  Therefore,  it  is  especially  important  in  this 
situation  to  calculate  the  value  of  e(u(P(8))}  for  the  solution  8 
obtained  in  order  to  verify  that  this  value  actually  is  positive  and 


reasonably  large.  It  also  would  be  advisable  to  rerun  the  procedure 


several  times  with  more  conservative  values  of  the  a ~  and  of  the 

J 

parameters  of  the  utility  function  to  see  if  less  risky  solutions 


yielding  a  higher  expected  utiJity  are  obtained. 


—  For  example,  if  N  independent  and  identically  distributed  invest¬ 
ments  are  approved,  then  cr(8)  =  -Jlicr.  whereas  p(8)  =  Np.. 

0  0 


Since  the  approximate  procedure  is  based  primarily  on  using  linear 
programing,  it  can  handle  a  large  number  of  constraints  with  ease, 
and  is  cor  utationally  feasible  far  problems  involving  an  extremely 
large  number  of  investment  proposals.  It  also  is  quite  well-suited 
for  conducting  an  extensive  sensitivity  analysis  (by  referring  to  the 
final  fractional  solution  rather  than  the  final  integer  solution). 
Another  advantage  is  that  it  can  be  readily  adapted  to  the  case  where 
same  of  the  investment  decisions  are  of  the  "how  much"  rather  than 
"yes-or-no"  type. 

The  exact  solution  procedure  does  not  enjoy  any  of  these  advantages 
of  the  approximate  procedure.  However,  it  does,  of  course,  have  the 
crucial  advantage  of  guaranteeing  an  optimal  solution  (with  respect 
to  the  objective  function  used)  for  any  problem  small  enough  for  it 
to  solve.  It  is  quite  efficient  for  problems  involving  approximately 
ten  investment  proposals  or  less.  However,  the  computation  time  grows 
rapidly  with  the  size  of  the  problem,  so  it  is  almost  certain  to  be 
computationally  feasible  only  up  to  approximately  fifteen  proposals. 

It  can  also  be  tried  with  at  least  a  small  hope  of  success  when  there 
are  say  20  to  25  proposals.  It  would  seem  that  the  presence  of 
mutually  exclusive  sets  of  investment  proposals  also  should  increase 
the  efficiency,  although  this  has  not  been  tested. 

When  it  is  assumed  that  P(B)  has  a  normal  distribution  and  that 
U(p)  follows  the  high  risk  aversion  model,  then  the  corresponding 
exact  expression  for  E{u(P(6) ) )  should  be  used  with  the  exact  pro¬ 
cedure.  For  other  .'Uses,  an  approximate  expression  based  on  a  Taylor 
series  expansion  of  E{u(P(B))}  is  available.  As  discussed  at  the 


62 


end  of  the  preceding  section,  this  approximate  expression  performs 
quite  well  (still  leading  to  the  optimal  solution  in  24  of  the  30 
test  problems  run  with  the  nigh  risk  aversion  model),  but  it  is 
sufficiently  imprecise  that  its  effect  should  not  be  ignored.  With 
the  high  risk  aversion  model,  it  usually  gives  somewhat  less  consider¬ 
ation  to  risk  than  it  should.  Therefore,  even  when  the  distribution 
of  P(6)  is  far  from  normal,  the  expression  that  is  exact  for  the 
normal  distribution  usually  should  be  used  instead  in  order  to  better 
approximate  the  appropriate  risk  tradeoffs.  The  only  exceptions  might 
be  distributions  that  are  skewed  substantially  to  the  rig.it,  since 
then  the  risk  for  a  given  p(6)  and  cr(6)  tends  to  be  somewhat  less 
than  with  the  normal  distribution.  When  the  basic  model  for  U(p)  is 
being  used,  the  approximate  expression  has  the  opposite  effect,  giving 
somewhat  more  weight  to  risk  than  is  warranted.  Since  no  exact 
expression  is  new  available  for  the  model,  it  is  then  especially 
important  to  conduct  sensitivity  analysis  to  check  for  alternative 
solutions  that  may  give  a  better  balance  between  risk  and  the  expected 
return  in  present  value. 

The  choice  between  the  two  alternative  models  for  U(p)  should 
be  based  both  on  their  inherent  properties  (as  discussed  in  Section  3) 
and  the  behavior  of  the  available  solution  procedures  with  them.  It 
appears  that  the  two  models  should  be  quite  comparable  in  their  evalu¬ 
ation  of  combinations  of  investments  having  only  a  moderate  degree  of 
risk,  but  the  high  risk  aversion  model  should  be  much  less  tolerant 
of  very  risky  combinations  than  the  basic  model.  In  effect,  the  high 
risk  aversion  model  has  a  threshold  level  of  tolerable  risk  beyond 


63 


which  it  increases  the  penalty  for  risk  very  rapidly  (due  to  the 
exponential  form  of  the  function).  Thus*  the  optimal  solutions  for 
the  two  models  (with  respect  to  exact  expressions  for  E(u(P(S)))) 
should  tend  to  be  fairly  comparable  (and  often  identical) ,  but  somewhat 
riskier  solutions  may  be  still  rated  good  by  the  basic  model  but  very 
poor  by  the  high  risk  aversion  model.  Because  of  the  conservative 
effect  of  using  the  approximate  expression  for  E(U(P(S))}  with  the 
basic  model,  the  solution  procedures  tend  to  occasionally  yield  somewhat 
less  risky  solutions  for  this  model  than  the  high  risk  aversion  model. 
This  effect  should  be  taken  into  account  in  choosing  between  the  two 
models.  However,  it  should  also  be  noted  that  this  tendency  may  be 
less  pronounced  (and  perhaps  even  reversed)  when  the  utility  function 
parameters  are  such  that  k  is  somewhat  less  than  one  (which  would 
be  common)  in  the  high  risk  {.version  model  rather  than  the  value  of 
one  used  for  almost  all  of  the  computational  experimentation.  A  final 
factor  in  making  the  choice  is  that  the  approximate  solution  procedure 
appears  to  be  at  least  slightly  more  effective  in  obtaining  good 
solutions  with  the  basic  model  than  with  the  high  risk  aversion  model. 

In  conclusion,  it  should  be  reemphasized  that  the  main  practical 
purpose  of  this  approach  to  analyzing  risky  interrelated  investments 
is  to  identify  and  quantitatively  describe  the  few  best  combinations 
of  investments  for  careful  evaluation  by  management.  Thus,  rather 
than  using  one  set  of  parameter  values,  one  utility  function,  and 
one  solution  procedure  to  find  "the  optimal  solution,"  it  usually  is 
important  to  conduct  t  sensitivity  analysis  on  the  data  to  identify 
those  solutions  that  are  best  under  some  plausible  set  of  estimates. 
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This  may  well  involve  using  more  than  one  utility  function  model  and 
solution  procedure ,  much  like  in  Tables  VI  and  VII.  This  process 
will  identify  a  small  number  of  alternatives  requiring  further  con¬ 
sideration,  and  will  provide  relevant  quantitative  measures  of 
performance  (u(&),  o(&),  the  resulting  risk  profile,  and  E{U(P(8)) }} 
for  these  alternatives  under  different  assumptions.  Management  can 
then  analyze  this  information,  as  well  as  other  relevant  factors  not 
incorporated  into  the  model,  and  make  its  decision  on  the  best  combi¬ 
nation  of  investment  proposals  to  approve. 
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APPENDIX 

Listing  of  Test  Problems 


This  appendix  presents  the  ten  basic  test  problems  that  were 
randomly  generated  as  described  by  Table  I  of  Section  4.  The  numbering 
of  the  problems  coincides  with  that  used  in  the  various  tables.  For 
each  problem,  the  diagonal  elements  of  "MU-MATRIX"  are  the  p.  in 
order,  and  the  off-diagonal  elements  are  the  p.,  .  The  diagonal 

JK 

elements  of  the  second  matrix  in  each  listing  are  the  s . ,  and  the 

'  J 

off-diagonal  elements  are  the  p...  The  third  matrix  lists  the  a 

JK  JK. 

and,  on  the  far  right,  the  b.. 
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